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Computer Assisted 
Methods for Diagnosing Diseases 



10 Cross Reference to Related Applications 

The present application is a continuation-in-part of co-pending U.S. 
Application Serial Number 08/315,851 filed on September 30, 1994, which 
is a continuation of U.S. Application Serial Number 07/990,772 filed on 
December 14, 1992, now abandoned, which is a continuation-in-part of 
15 U.S. Application Serial Number 07/964,486 filed on October 21, 1992, 

now abandoned, which is a continuation of U.S. Application Serial Number 
07/806,980 filed on December 12, 1991, now abandoned. 



20 Technical Field 

The present invention relates to methods for diagnosing, screening 
or prognosing diseases. More particularly, the present invention relates to a 
method for diagnosing, screening or prognosing diseases in humans or 
animals, and for determining the severity and cause of the disease, by 

25 measuring blood serum levels of specific, predetermined blood constituents 

and then calculating a diagnostic index based on the relationships between 
those blood constituents. Optionally, demographic information, such as the 
age, race or sex of the patient, as well as the patient's medical history can be 
factored into the diagnosis. 

30 The present invention further relates to an in vitro, serum -based 

computer assisted method for diagnosing, screening or prognosing 
diseases, utilizing a neural network to obtain a diagnostic index. In 
preferred embodiments of the present invention, the method is used to 
diagnose osteopenia and certain cancers, including but not limited to 

35 ovarian, breast, testicular, colon and prostate cancer. 
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Background of the Invention 

As used herein, the term "disease" is defined as a deviation from the 
normal structure or function of any part, organ or system of the body (or 
any combination thereof). A specific disease is manifested by characteristic 
symptoms and signs, including both chemical and physical changes. 
Certain characteristic signs and symptoms can be quantitated by clinical 
chemical analysis to yield important diagnostic information. For purposes 
of this application, the quantifiable signs, symptoms and/or analytes in 
biological fluids characteristic of a particular disease are defined as 
"biomarkers" for the disease. Current diagnostic methods depend on the 
identification and evaluation of these biomarkers, both individually and as 
they relate to one another. 

Generally, the pathological process involves gradual changes that 
become apparent only when overt change has occurred. In many instances, 
pathological changes involve subtle alterations in multiple biomarkers. It is 
uncommon that a single biomarker will be indicative of the presence or 
absence of a disease. It is the pattern of those biomarkers relative to one 
another and relative to a normal reference range, that is indicative of the 

presence of a disease. 

When individual biomarkers do not show a predictable change and 
the patterns and interrelationships among the biomarkers viewed collectively 
are not clear, the accuracy of a physician's diagnosis is significantly 
reduced. Also, as the number of biomarkers significant for a particular 
disease increases, the number of diagnostic patterns increases, 
correspondingly, pattern recognition and accuracy of a physician's 

diagnosis decreases. 

For example, the following biomarkers collectively show 
characteristic changes in the presence of osteoporosis: calcium, phosphate, 
estradiol (follicular, mid-cycle, luteal, or post-menopausal), progesterone 
(follicular, mid-cycle, luteal, mid-luteal, oral contraceptive, or over 60 
years), alkaline phosphatase, percent liver-ALP, and total intestinal- ALP. 
After measuring these biomarkers, a diagnosing physician would next 
compare the measuring to a normal reference range. While some of the 
biomarkers may fall outside the normal reference range, others may fall 
clearly within the normal reference range. In some circumstances, all of the 
biomarker values may fall within a normal reference range. Presented with 
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such data, a physician may suspect that a patient has undergone some bone 
loss, but will be unable to reach a conclusive and meaningful diagnosis as to 
the presence of the disease osteoporosis. 

The characteristic changes in biomarkers associated with some 

5 diseases are well documented; however, the quantitative interpretation of 

each particular biomarker in diagnosing a disease and determining a 
prognosis is not well established. The difficulties inherent in formulating a 
diagnosis from the analysis of a set of laboratory data is best illustrated by 
looking closer at conventional diagnostic methods for a specific disease. A 

10 discussion of the disease osteoporosis follows. 

The term "osteopenia" as used herein means any decrease in bone 
mass below the normal. The term "osteoporosis" as used herein means a 
specific form of generalized osteopenia characterized by a decrease in bone 
density, low bone mass, and microarchitectural deterioration of bone tissue. 

15 Osteopenia encompasses a group of diseases with diverse etiologies 

typified by reduction in bone mass per unit volume to a level below that 
which is necessary for adequate mechanical support. Osteoporosis is the 
result of the gradual depletion of the inorganic portion of the skeleton and 
can be caused by any number of factors. Primary osteoporosis is an age 

20 related disorder that is particularly common in women and is characterized 

by decreased bone mass in the absence of other recognizable causes. 
However, osteoporosis occurs in both men and women. In women it is 
recognized usually at the 5th or 6th decade, following menopause. In men 
osteoporosis is often recognized around their 6th or 7* decade. 

25 The following is a partial list of some of the categories of individuals 

at risk for developing osteoporosis: 

Postmenopausal women 

Cigarette smokers 
30 Heavy users of alcohol 

Users of a variety of drugs, such as steroids 

Female runners and ballet dancers 

Male marathoners consuming two few calories 

Bulimics and anorexics 
35 People with poor diets 

People allergic to dairy products 
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People affected with cancer 

Fair and slim women 

All men and women over the age of 65. 

In addition to being female, the three most significant risk factors are 
poor diet, lack of exercise, and being postmenopausal. Other risk factors 
which are associated with osteoporosis include Caucasian or Oriental 
ancestry, a fair complexions, and a family history of osteoporosis. 

The onset of osteoporosis may be insidious or sudden, following 
trauma. The most common complaint associated with osteoporosis is back 
pain. Eventually, the pain may spread to the pelvis, the thorax, and the 
shoulders. In the spine, the vertebrae can compress, and the back can take 
on a "bent" appearance. Conditions such as kyphosis (humpback) or 
scoliosis may occur. If the spine becomes deformed, other body parts can 
be affected as well. For example, the ribs can be pushed against the pelvis, 
or the stomach can be pushed into the pelvis. In addition to spinal 
problems, osteoporosis can also lead to fractures of the hip, wrist, and nhs 
These fractures can occur with only slight trauma and sometimes with no 
trauma at all. Mazess B., et aL 9 "Bone Density of the Radius, Spine, and 
Proximal Femur in Osteoporosis," J. of Bone and Mineral Research, Vol 
3, pgs. 13-18, (1988); Riggs B. L., et aL 9 "Involutional Osteoporosis". 
New Engl. J. Med., Vol. 314, pgs. 1676*1686, (1986). The changes 
associated with osteoporosis are gradual so osteoporosis is often not 
detected in its early stages. 

From the research perspective, osteoporosis can be classified as 
either primary or secondary to another disease. Primary osteoporosis is 
further classified as juvenile, idiopathic, postmenopausal (Type I) and 
involutional (Type II). It is now understood that accelerated bone loss 
occurs with cessation of menstruation at the time of menopause and in 
women who have amenorrhea as a result of prolactin producing pituitary 
tumor, anorexia nervosa, or intense long-distance running associated with 
undernourishment. These situations are all accompanied by estrogen 
deficiency which is likely to be a major determinant of the accelerated bone 
loss. Bone loss also occurs when estrogen therapy is withdrawn. 

Bone consists primarily of an extracellular matrix containing (by 
weight) approximately 35% organic and 65% inorganic components. The 
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cells of bone represent a minor component of the bone constituents, yet they 
carry out a major portion of the function of the skeletal system. Bone cells 
help maintain serum calcium concentration within a narrow range to regulate 
mineral homeostasis while also being responsible for the continuous 
5 formation and resorption of the extracellular matrix, allowing response of 

the skeletal system to the mechanical forces resulting from physical activity. 

The major organic component of the extracellular matrix is collagen. 
This protein has a rigid rodlike structure and is composed of three alpha 

10 chains held together in helical fashion by covalent and noncovalent forces. 

Multiple collagen molecules form fibrils, and these fibrils in turn are 
arranged in bundles or fibers. It is the bundles or fibers of collagen that can 
be seen in the light microscope as layers or linear arrays. 

Non-collagen components of bone compose a very small portion of 

15 the organic matrix of the skeleton. This minor fraction consists of proteins, 

glycoproteins, lipids and mucopolysaccharides. Only a few of these 
components have been carefully identified and characterized. Two of the 
proteins that have been isolated and studied are osteonectin and osteocalcin. 
Calcium and phosphorus are the main components of the inorganic 

20 portion of the skeleton. Initially, calcium and phosphorus are deposited as 

amorphous salts but later undergo rearrangements into a crystalline structure 
that resembles hydroxyapatite (Caio(P04>6(OH)2). Several other ions, 
including Na, K, Mg, and CO3 in varying proportions, may be found in the 
skeletal hydroxyapatite. If there has been fluoride intake, there will also be 

25 fluoride in the hydroxyapatite. 

Although considerable effort has been expended in studying the 
mechanism of bone mineralization, it is not fully understood. Several 
theories have been suggested to explain the available data. One theory is 
that calcium and phosphorus ions are present in the extracellular fluid in 

30 amounts exceeding the solubility product of [Ca] x [P], These ions are 

kept from precipitating by inhibitors of calcification such as pyrophosphate. 
Because osteoblasts contain large amounts of alkaline phosphatase, it has 
been speculated that the activity of this enzyme facilitates mineralization by 
cleaving the phosphate groups, thus altering the Ca:P ratio in the sites of 

35 calcification. 
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In the clinical syndrome of osteoporosis, the reduction in bone mass 
can be attributed to osteopenia because of a dietary deficiency or absorption 
interference of proteins or Vitamin C. It can also be the result of a deficient 
stress stimulus. Osteopenia can also be caused by osteomalacia, a failure of 
proper mineralization of osteoid resulting from bone calcium or phosphorus 
deficiency or both. It can also be caused by insufficient absorption from the 
intestine due to lack of calcium or a resistance to the action of Vitamin D due 
to failure of its conversion to the biologically active forms, 25 
hydroxychole-calciferol and 1.25 dihydroxychole-calciferol, formed by the 
liver and kidney, respectively. In addition, it may well be caused by an 
abnormal rate of osteolysis due to parathyroid hormone stimulation of 
osteoclastic activity (osteocasts are hematopoitic in origin, arising from the 
migration of monocytes to bone). In reality, most cases of osteoporosis, 
when carefully analyzed, reveal evidence of the causes of osteopenia. 

Chemical analysis of blood may reveal calcium, phosphorus, and 
alkaline phosphatase within the normal range. However, an isoenzyme of 
alkaline phosphatase may be significantly increased. Increased bone 
resorption seen in osteoporotic patients, which occurs as a result of the 
action of osteoclasts, usually involves the dissolution of both minerals and 
organic matrix eventually leading to increased excretion of urinary 
hydroxyproline. Serum estradiol which is secreted almost entirely by the 
ovary is significantly decreased in these patients. This observation is further 
corroborated where it has been demonstrated that exogenous estrogen 
therapy in perimenopausal women does delay the onset of post menopausal 
osteopenia. Most experts agree that estrogen appears to decrease bone 
resorption. Weiss N. S., et aU "Decreased Risk of Fractures of the Hip 
and Lower Forearm with Postmenopausal Use of Estrogen," N. Engl. J. 
Med., Vol. 303, pgs. 1195-1198 (1980); Ettinger et al, "Long Term 
Estrogen Replacement Therapy Prevents Bone Loss and Fractures," Ann. 
Intern. Med., Vol. 102, pgs. 319-324, (1985). 

It is believed that blood levels of calcium are maintained without 
increasing calcium loss from the bone as long as there are normal amounts 
of estrogen and parathyroid hormone present. It is reasonably excepted that 
estrogen antagonizes the effect of parathyroid hormone. Estrogen 
deficiency, as seen in peri- and postmenopausal women, results in an 
increase in the sensitivity of bone to parathyroid hormone. This antagonistic 



WO 96/12187 



PCT/US95/01379 



7 



relationship eventually leads to an increase in the resorption of the bone and 
contributes to the development of Type I Osteoporosis. Type II 
Osteoporosis is characterized by reduced calcium absorption, which in turn 
results in increased secretion of parathyroid hormone. Type n Osteoporosis 

5 occurs in older individuals and is associated with wedge fractures of the 

spine and fractures of the hip. A theory recently advanced by Barnhill, et al 
suggests that evidence of decreased serum estradiol and increased 
lymphocyte alkaline phosphatase, represents an activated immune system in 
osteopenic postmenopausal women. This finding suggests that "uncovered 

10 estrogen receptors" may induce an immune reaction which is responsible for 

one form of osteoporosis. Barnhill S.D., et al, "Osteoporosis: A Possible 
Autoimmune Etiology," Ann. of Clin. Lab. Sci., Vol. 17, pgs. 255-256, 
(1987). 

In 1987, Barnhill, et al proposed an autoimmune etiology for some 

15 forms of osteopenia. In that publication, Barnhill et al showed that 

lymphocyte-derived alkaline phosphatase was present in the blood of 90% 
of severely osteopenic women. See Barnhill etal, "Osteoporosis: A 
possible Autoimmune Etiology", Ann. of Clin. Lab. Scl, Vol. 17, pgs. 
255-256, (1987)). The concept of a lymphocyte-derived alkaline 

20 phosphatase is further described by Griffiths, et al See Griffiths, J., et al, 

"Separation and Identification of Alkaline Phosphatase Isoenzymes and 
Isoforms in Serum of Healthy Persons by Isoelectric Focusing", Clin 
Chem, Vol. 32, pgs. 2171-2177, (1987). Thus, osteoporosis is well suited 
to diagnosis by interpretation of laboratory data. 

25 Osteoporosis poses a major health problem in the United States, not 

only for those persons who are already affected but for individuals whose 
diet, life style, and body build make them more likely to develop 
osteoporosis as they age. Postmenopausal osteoporosis is a common 
disorder that results in substantial morbidity and mortality. As many as 25 

30 million American women suffer from osteoporosis. Osteoporosis is 

responsible for approximately 250,000 femoral neck fractures annually in 
the United States, and these fractures are associated with a 20% mortality 
rate within 6 months. The risk is particularly high among the elderly, who 
also tend to lose bone as a result of the aging process. Wu K., et al, "Bone 

35 Resorption Rates in Physiological, Senile, and Postmenopausal 

Osteoporosis." J. Lab Clin. Med., Vol. 69, pg. 810, (1967). 
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It is estimated from a survey of medical clinics that of those 
individuals living to age 85, 32% of the women and 17% of men will 
fracture hips, weakened by osteoporosis. In addition to the pain and 
suffering caused by these fractures, the monetary cost is great, accounting 
for well over 3.8 billion dollars per year in treatment for fractures of 
osteoporosis. Moreover, six to eight months following hip fracture, about 
50% of the osteoporotic patients are in need of assistance with activities of 
daily living and about 25% require nursing home care. Only about 25% of 
the patients fully recover. In view of the costs associated with osteoporosis, 
as measured in dollars and in human suffering, osteoporosis has 
increasingly been perceived as a serious and disabling disease, warranting 
substantial involvement on the part of clinical investigators, governmental 
agencies, and pharmaceutical industries to develop and evaluate potential 
treatments and early detection techniques. 

The National Institutes of Health 1984 Consensus Conference 
stimulated those interested in bone density with its recommendations for 
"defining persons at risk, and developing safe, effective and low cost 
strategies for fracture protection." See, Office of Medical Applications of 
Research National Institutes of Health, "Osteoporosis: consensus 
conference," JAMA, Vol. 252, pgs. 799-802, (1984). Since 1984, there 
has been a dramatic change in the way physicians view osteopenia. With 
the development and popularization of highly sensitive techniques such as 
DPA and quantitative computed tomography, physicians are now capable of 
measuring the density of the proximal femur and the lumbar vertebrae. (For 
a general review of the methods currently available to measure bone density, 
see Avioli, L.V., ed., "Metabolic Bone Disease and Clinically Related 
Disorders", W.B. Saunders Company (1990)). 

An early decrease in bone mass can be measured by non-invasive 
assessment of the skeleton by four widely available methods, including 
single photon absorptometry, dual photon absorptometry, dual-energy x-ray 
absorptometry, and quantitative computed tomography. A brief discussion 
of these diagnostic methods follows. 

A device called a single-photon absorptometer (SPA) is used to 
measure bone mineral content, primarily in the forearm and wrist. The heel 
bone can also measured using SPA because the heel bone is thought to be a 
predictor of bone loss in the spine. SPA measures primarily cortical bone, 
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which is also affected by osteoporosis, though not to the same extent as 
trabecular bone. 

The technique of dual-photon absorptometry (DPA) provides a 
measurement of the total cortical and trabecular mineral content of the hip 

5 and spine. DPA uses less radiation than conventional x-rays; but, a scan of 

the spine using DPA still exposes the patient to approximately one tenth of 
the radiation that results from a routine chest x-ray. 

Dual-energy x-ray absorptometry (DXA) provides a measurement of 
the amount of bone tissue in the hip and spine. This technique is now used 

10 routinely because it is faster than DPA. 

Unlike DPA and DXA, quantitative computed tomography, more 
commonly called a CAT scan, can measure the density of either bone or just 
the trabecular portion. CAT scans unfortunately expose patients to higher 
doses of radiation than any of the other techniques. 

15 Radiographic absorptometry (RA) is a method for non-invasive 

measurement of bone mineral x-rays of the hand. Radiographs, taken with 
a standard x-ray machine, are sent to a central laboratory for computer- 
controlled analysis. 

As osteopenia progresses, both the number and thickness of the 

20 trabecular units decrease, causing fragility of bone and an increased risk of 

fractures. The radiographic manifestations of osteoporosis reflect the 
deficiency of the organic matrix and parallel the gross pathologic findings. 
The cortices are thin and the trabeculae fine and sparse; the skeletal 
structures are therefore more radiolucent than normal. The disease process 

25 eventually affects almost all of the skeletal structures. The principal areas of 

demineralization are the spine and pelvis, especially in the femoral neck and 
head. Demineralization is less marked in the skull and extremities. 
Although useful to detect breaks in the bone, ordinary x-rays are not 
sensitive enough to detect osteoporosis until a large amount of bone tissue 

30 has already been lost, generally from 25% to 40%. By the time 

osteoporosis can be identified by x-ray techniques, the condition is 
advanced. 

Current standard diagnostic techniques, are not effective for early 
detection of osteoporosis. Changes seen in osteoporosis are very gradual, 
35 and often go undetected in the early stages of the disease. Osteoporosis is 

often not detected in its early stages because bone mass must be decreased 
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by about 30% to 40% before it is apparent using standard x-ray diagnostic 
techniques. Preventing osteoporosis by detecting early bone loss is far 
better than identifying the disease at relatively advanced stages and 
subsequently attempting to prevent its progression. Once major 
deterioration has occurred and gaps exist between the ends of fractured 
trabecular, no current treatment can be expected to restore the lost bone. 
Thus, therapeutic efforts must be directed toward prevention and early 
recognition of the progressive disease so treatment can be instituted before 
essentially irreversible structural damage ensues. Cummings 5./?., et al, 
"Should Perimenopausal Women Be Screened for Osteoporosis?", Ann. 
Int. Med., Vol. 104, pgs. 745-751, (1986); Courpron P., "Bone Tissue 
Mechanisms Underlying Osteoporosis," Orthop. Clin. North Am., Vol. 
12, pg. 513, (1981); Frost H. M., "Mechanical Determinants of Bone 
Modeling," 7. Metabol Bone. Dis. Rel Res., Vol. 4, pg. 217, (1983). 

One of the problems with the current methods for diagnosing 
osteoporosis is that the procedures do not give any information about the 
underlying cause of the osteoporosis, making it difficult to prescribe an 
appropriate course of treatment for the patient. For example, a common 
cause of postmenopausal osteoporosis is an estrogen deficit, which x-ray 
techniques cannot measure. Another problem inherent in the current 
diagnostic methods for osteopenia is that all of the current methods require 
expensive, sophisticated medical instrumentation to perform the bone 
density measurements. Additionally, patients must be exposed to x-rays. 
This makes a general screening of high risk populations impractical due to 
the expense and unavailability of the necessary instrumentation to the 
average clinic. 

In view of the difficulties associated with extracting a diagnosis from 
the laboratory data for a set of predictive biomarkers, there is need for 
automated diagnostic systems that are capable of complex pattern 
recognition. There have been several attempts at using computational 
models to achieve pattern recognition in diagnostics. One of the most 
popular computational method for making diagnoses from multivariate 
laboratory data has been discriminate function analysis. However, 
diagnostic systems that rely exclusively on classical pattern recognition 
technology (geometric, syntactic, template, statistical) are not effective for 
evaluating the characteristic biomarker patterns of many disease states. 
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There is no clear set of rules that accurately describes how to analyze a set of 
biomarkers to reach a diagnosis. 

In recent years, neural networks have been gaining popularity as a 
means for recognizing and analyzing subtle diagnostic patterns in 

5 multivariate laboratory data. Neural networks possess the ability to discern 

patterns and trends too subtle or too complex for humans and conventional 
computational methods to identify. While humans can't easily assimilate 
more than two or three variables at once, neural networks can perceive 
correlations among hundreds of variables. Examples of areas in which 

10 neural networks have been explored for their value in clinical diagnosis 

and/or prognosis include: 

• psychiatry (See Mulsant, B.H., "A Neural Network as an 
Approach to Clinical Diagnosis", MD Computing, Vol. 7, pp. 25-36 
(1990)); 

15 ♦ autism (See Cohen, I., et al, "Diagnosing Autism: A Neural 

Net-Based Tool", PCAI, pp. 22-25 (May/June 1994); 

• pediatric radiology (See Boone, L.M., et al., "Neural Networks 
in Radiologic Diagnosis. I. Introduction and Illustration", Invest Radiol, 
Vol. 25, pp. 1012-1016, (1990) and Gross, G.W., et al, "Neural 

20 Networks in Radiologic Diagnosis. II. Interpretation of Neonatal Chest 

Radiographs", Invest Radiol, Vol. 25, pp. 1017-1023 (1990)); 

• breast cancer (See Astion, M.L., et al, "Application of Neural 
Networks to the Interpretation of Laboratory Data in Cancer Diagnosis", 
Clin. Chem., Vol. 38, No. 1, pp. 34-38 (1992); Yuzheng, W., et al, 

25 "Artificial Neural Networks in Mammography: Application to Decision 

Making in the Diagnosis of Breast Cancer", Radiology, Vol. 82, pp. 81-87 
(1993); Kappen, H.J., et al, "Neural Network Analysis to Predict 
Treatment Outcome", Annals of Oncology, Vol. 4, Supp. 4, pp. S31-S34 
(1993); and, Ravdin, P.M., et al, " A practical application of neural 

30 network analysis for predicting outcome of individual breast cancer 

patients", Breast Cancer Research and Treatment, Vol. 22, pp. 285-293 
(1992)); 

• ovarian cancer (See Wilding, P., et al, "Application of 
backpropogation neural networks to diagnosis of breast and ovarian 

35 cancer", Cancer Letters, Vol. 77, pp. 145-153 (1994)). 
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• thyroid disease (See, Sharpe, P.K. t et. ah\ "Artifical Neural 
Networks in Diagnosis of Thyroid Function from in Vitro Laboratory 
Tests," Clin. Chem., Vol. 39, No. 11, pps. 2248-2253 (1993)). 

• cervical cancer (See U. S. Patent No. 4,965,725 to Rutenberg); 

5 and, 

• cardiology (See U.S. Patent No. 5,280,792 to Leong et al and 
Furlong, J.W., "Neural Network of Serial Cardiac Enzyme Data: A Clinical 
Application of Artifical Machine Intelligence", Clin. Chem., Vol. 96, No. 1, 
pp. 134-141 (July 1991). 

10 Neural networks are capable of pattern recognition particularly suited 

to making diagnoses. Unlike current methods for arriving at a diagnosis 
from a logical set of rules, neural networks do not require explicit encoding 
of process knowledge in a set of rules. Neural networks learn from 
examples. 

15 The present invention provides a new method for diagnosing 

diseases and screening patients at risk for a disease, utilizing laboratory data 
to obtain a diagnostic index. It further provides a method for early detection 
and diagnosis of a disease. As applied to osteoporosis, it generates a new 
level of interest in screening individuals who are at risk for osteoporosis 

20 because this novel diagnostic procedure provides information about the 

underlying cause of the osteopenia. The present invention also provides a 
simple, inexpensive, and rapid method of estimating bone density using 
serum blood levels of specific, predetermined blood constituents. As 
applied to cancer, it provides a simple, inexpensive, accurate and rapid 

25 method for detecting cancer, resulting in earlier and less invasive treatment. 

Summary of the Invention 

The present invention is directed to an in vitro, biological fluid-based 
method for diagnosing, screening or prognosing for diseases in which 

30 laboratory data are manipulated to obtain a diagnostic index. The term 

"biological fluid" includes, but is not limited to, blood, urine, saliva, 
cerebral spinal fluid, synovial fluid, and tears. In the preferred 
embodiments of the present invention, the method of the present invention is 
used to provide a of cancer and osteoporosis. 

35 The present invention is a method for diagnosing, screening or 

prognosing a disease in a human or animal comprising the steps of 
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measuring the concentrations of a predetermined set of biomarkers known to 
be associated with the disease from a biological fluid from the human or 
animal; scaling the digitized values of the analytes; and sending the scaled 
values to a trained neural network, whereby the output values from the 

5 neural network tend toward the upper value when the human or animal has 

the disease and the output values tend toward the lower value when the 
human or animal does not have the disease. 

The present invention also comprises an apparatus for diagnosing, 
screening or prognosing a disease in a human or animal comprising a means 

10 for digitizing the concentrations of a predetermined set of biomarkers known 

to be associated with the disease from a biological fluid from the human or 
animal; a means for scaling the digitized values; and a trained neural network 
coupled to the digitizing and scaling means for generating network output 
values between an upper and lower value, whereby the output values tend 

15 toward the upper value when the human or animal has the disease and the 

output values tend toward the lower value when the human or animal does 
not have the disease. 

The present invention also comprises a method for diagnosing, 
screening or prognosing a disease in a human or animal comprising the 

20 steps of measuring the concentrations of a predetermined set of biomarkers 

known to be associated with the disease from a biological fluid from the 
human or animal, scaling the digitized values of the analytes, and 
introducing the scaled values to a first trained neural network, whereby the 
output values from the first neural network tend toward the upper value 

25 when the human or animal has the disease and the output values tend toward 

the lower value when the human or animal does not have the disease; and 
sending the output value from the first neural network and a second set of 
predetermined biomarkers, which could include one or more of the 
biomarkers in the first set of predetermined biomarkers, to a second trained 

30 neural network, whereby the output values from the second neural network 

tend toward the upper value when the human or animal has the disease and 
the output values from the second neural network tend toward the lower 
value when the human or animal does not have the disease. 

In accordance with the first embodiment of the present invention, a 

35 trained neural network is utilized to determine a diagnostic index 

corresponding to the presence and severity of a disease by analyzing a set of 
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predetermined biomarkers for that disease. In accordance with the 
invention, the concentrations, or in some cases, the presence of certain 
biomarkers related to the incidence of a particular disease are determined for 
a patient These values are then scaled and the scaled values are then sent to 
a trained neural network to yield a diagnostic index. A neural network is 
trained by introducing a population of patients in which a disease state is 
known, along with the biomarker values for those patients and "teaching" 
the neural network to recognize the patterns in the biomarkers. After the 
neural network is trained, biomarker values from patients with unknown 
disease states are introduced to the trained neural network. The neural 
network then processes the information to produce a value corresponding to 
a diagnosis of the presence or absence and the severity of a particular 
disease. 

The inventors propose that the artificial neural network, especially 
the multi-layer feedforward network, may, through their weighi 
connections, correspond to biomarker patterns that are important for 
categorizing diseases. Additionally, the inventors propose that the neural 
network is able to identify unique patterns of biomarkers associated with a 
variety of disorders that may help to classify borderline cases that do not 
appear to fit into either a malignant or benign pattern. 

A second embodiment of the present invention involves a two sicp 
analysis of the biomarkers by neural network. This avoids the bias created 
by a dominant predictive variable when training a network. The dominant 
biomarker or predictive variable is excluded from the first analysis by neural 
network and then included in a second analysis by neural network. For 
example, if age is thought to be the dominant predictive variable in the 
diagnosis of osteoporosis, that variable is not included in the training of the 
first neural network, and the training data set is limited to the other selected 
biomarkers. After obtaining a diagnostic index using the first set of 
biomarkers, a second neural network is trained using the diagnostic index 
and the entire set of input variables, including age, to yield another 
diagnostic index. The final diagnostic index is a composition of an artificial 
neural network generated index and results from heuristic analysis using 
other non-numerical patient information. 

The present invention further comprises an apparatus capable of 
diagnosing, screening or prognosing a disease comprising a sample 



WO 96/12187 



PCT/US95/01379 



15 



receiving means, a sample detecting means capable of detecting the quantity 
of one or more biomarkers in a biological fluid, and, after digitalizing the 
quantity of analyte in the biological fluid, scaling the digitized values and 
introducing the scaled data into a neural network which yields a therapeutic 

5 index on a printer or, optionally, on a video display. 

A third embodiment of the present invention includes a method for 
diagnosing and determining the severity and underlying cause of osteopenia 
using blood concentrations of a specific, predetermined set of blood 
constituents. The method comprises the steps of (1) measuring the severity 

10 of the disease in a set of humans or animals with varying severity of disease 

by a standard method, (2) assigning the severity of disease a numerical 
value on a severity scale, the scale being from no disease to severe disease. 
(3) measuring the blood concentrations of the predetermined set of blood 
constituents in the set of humans or animals with varying severity of the 

15 disease, and then (4) calculating a numerical relationship between the set of 

blood concentrations and the severity of disease. 

The present invention also comprises a simple and rapid method of 
determining severity of osteopenia in a patient In a preferred embodiment 
the method comprises determining the serum level of the following scrum 

20 constituents: calcium, phosphate, total alkaline phosphatase, an alkaline 

phosphatase isoenzyme, estradiol, and progesterone. The alkaline 
phosphatase isoenzyme is preferably t-lymphocyte derived alkaline 
phosphatase or blood, liver or intestinal alkaline phosphatase isoenzyme . 
The results of these tests are then introduced into an algorithm. Optionally, 

25 the age of the patient may also be factored into the equation. The bone 

density coefficient that is calculated by the algorithm correlates to a very 
high degree to bone density as measured by standard methods, such as 
radiographic absorptometry, quantitative computed tomography, dual 
photon absorptometry and direct measurement of bone density. The bone 

30 density coefficient that is measured is then compared to an osteopenic 

severity scale. 

Using the six serum constituent concentrations, the present invention 
can be used to determine the osteopenic state of a patient as well as give an 
indication of the underlying cause of the osteopenia. The present invention 
35 can be correlated to any method of measuring bone density simply by 
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recalculating the coefficients in the algorithm using a multiple linear 

regression analysis. 

Accordingly, it is an object of the present invention is to provide a 
method for diagnosing, screening or prognosing and determining the 
severity of a disease using the concentrations and values of a predetermined 

set of biomarkers. 

It is a further object of the present invention to provide a more 
accurate method for screening, prognosing and diagnosing breast cancer 
using biological fluid-based laboratory data. 

It is another object of the present invention to provide a more 
accurate method for screening and diagnosing ovarian cancer using 
biological fluid-based laboratory data. 

Yet another object of the present invention is to provide a more 
accurate method for screening, prognosing and diagnosing colon cancer 
using biological fluid-based laboratory data. 

Furthermore, it is an object of the present invention to provide a 
more accurate method for screening, prognosing and diagnosing prostate 
cancer using biological fluid-based laboratory data. 

Additionally, it is an object of the present invention to provide a 
more accurate method for screening, prognosing and diagnosing testicular 
cancer using biological fluid-based laboratory data. 

Another object of the present invention is to provide a method for 
diagnosing cancer which will provide a better understanding of the probable 

cause of the cancer. 

Another object of the present invention is to provide a diagnostic test 
for cancer which can be used to screen large numbers of individuals. 

Another object of the present invention is to provide a simple and 
rapid chemical test for diagnosing osteoporosis. 

Another object of the present invention is to provide a test for 
osteoporosis which will also give information as to the underlying cause of 
the osteopenic condition. 

Another object of the present invention is to provide a diagnostic test 
for osteoporosis which can be used to screen large numbers of individuals. 

Yet another object of the present invention is to provide a method for 
diagnosing osteoporosis and determining the underlying cause of the 
osteopenia without having to subject the patient to radiation. 
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These and other objects, features and advantages of the present 
invention will become apparent after a review of the following detailed 
description of the disclosed embodiment and the appended claims. 

5 Brief Description of the Figures 

Fig, 1 illustrates a feed forward neural network having multiple 

outputs. 

Fig. 2 illustrates a feed forward neural network having a single 

output. 

10 Fig. 3 is an equation illustrating the mathematical relationship 

between the input and output of a neuron. 

Fig. 4 is a schematic illustration of the second preferred 
embodiment of the present invention. 

Fig. 5 is a graph showing the change in the neural network 
15 diagnostic index obtained by the method of the present invention over time 

for a ovarian cancer patirent 

Fig. 6 is a graph showing the change in the neural network 
diagnostic index obtained by the method of the present invention over time 
for a prostate cancer patirent 
20 Fig. 7 is a graph showing the change in the neural network 

diagnostic index obtained by the method of the present invention over time 
for a colon cancer patirent 

Fig. 8 is a graph showing the change in the neural network 
diagnostic index obtained by the method of the present invention over time 
25 for a breast cancer patient. 

Fig. 9 shows the patient data used to develop and test the 
osteoporosis neural network diagnostic system. 

Fig. 10 shows the patient data used to test the osteoporosis neural 
network diagnostic system. 
30 Fig. 11 is a program list of the algorithm of the third embodiment 

of the present invention. 

Fig. 12 shows the discrepancies between the DP A measurements 
and the measurements obtained using the method of the present invention 

for individual subjects. 
35 Fig. 13 shows the correlation of the results obtained using the 

present invention and DPA. 
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Fig. 14 shows the normality of values predicted by the present 
invention. 

Fig. 15 shows the normality of residuals (differences between DP A 
measurements and values predicted according to the method of the present 
5 invention), as evidenced by the arrangement of points along a diagonal 

straight line. 

Fig. 16 shows that the variance of residuals. 
Fig. 17 shows the independence of residuals. 
Fig. 18 shows that virtually all observed DPA measurements are 
10 predicted by the method of the present invention. 

Fig. 19 shows the training data used to construct the prostate 
cancer neural network prognostic system. 

Fig. 20 shows the relationships between the input neurons, hidden 
layer and output neuron of a neural network in one embodiment of the 
15 present invention. 

Detailed Description 

The first embodiment of the present invention is directed to an in 
vitro, serum-based computer assisted method for screening, prognosing and 

20 diagnosing diseases utilizing a neural network to obtain a conclusive 

diagnosis. The present invention can be adapted to existing diagnostic 
devices that have a collection means, a sample detecting means capable of 
detecting the quantity of an analyte in a biological fluid and a means of either 
printing or displaying the results of the tests on video display means. 

25 The inventors have discovered that biomarkers collectively alter in 

response to a disease process, and collectively constitute a new diagnostic 
biomarker with better disease predictability than the individual biomarkers. 
When the biomarkers are processed and analyzed as a group to yield a single 
diagnostic index, the sensitivity and specificity of the diagnosis is increased, 

30 making it possible for a physician to detect the presence of a disease earlier 

and with greater precision, or estimate a prognosis with greater precision, 
than by analysis of the individual biomarkers. 

In accordance with the method of this invention, blood is first 
collected from a patient, and then centrifuged to separate the red blood cells 
35 from the serum. It should be understood that biological fluids other than 

blood can be used in practicing the present invention. Next, the serum is 



WO 96/12187 



PCT/US95/01379 



19 



analyzed, using standard laboratory techniques, to determine the 
concentrations, or in some cases the presence or absence, of specific 
predetermined biomarkers related to the incidence of a particular disease. It 
is to be understood that this process can be carried out automatically in 

5 conventional diagnostic machines. For purposes of illustration, a 

description of the methods for obtaining the values for the biomarkers for 
osteopenia is provided elsewhere in this section. 

The biomarkers relied upon to diagnose a disease by the method of 
the present invention must be predictive of the suspected disease and must 

10 be statistically significant for analysis by a neural network. The selection of 

biomarkers that offers statistically significant discriminating power in the 
diagnosis of disease involves several steps. First an inventory of 
biomarkers that have shown certain relevancy in the diagnosis of the disease 
of interest must be conducted. In general, only the biomarkers that reflect 

15 different aspects of the disease process or other diagnostic information need 

to be included. Second, the selected biomarkers need to have a reasonable 
diagnostic value in terms of sensitivity, specificity, and positive and 
negative predictive powers. The design and implementation of experimental 
protocol from which the biomarkers are developed and evaluated should 

20 also be considered. Third, if the number of candidate biomarkers is large, a 

formal discriminating power analysis may be conducted. However, many 
of the standard statistical analysis methods may not be adequate for highly 
nonlinear classification problems. If the number of candidates are not too 
large, they may be all included in the initial attempt of neural network 

25 training. If one or several of the input biomarkers to the network are 

irrelevant to the classification decision making process, it will be reflected in 
the network connection weights of the trained neural networks. These 
values may then be removed from the biomarker set for a particular disease. 
Other methods for evaluating the statistical significance of a biomarker 

30 selected for analysis by neural network and selecting biomarkers for training 

a neural network are well known in the art. 

Biomarkers which meet the criteria delineated above, namely, they 
are predictive of a particular disease and statistically significant for analysis 
by neural network, are identified below for a sample of diseases including 

35 ovarian cancer, prostate cancer, colon cancer, breast cancer, testicular cancer 

and osteoporosis. 
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Ovarian Cancer 


Prostate Cancer 1 


Colon Cancer 


LASA-P® 


LASA-P® 


LASA-P® 


CA125 


PAP 


CA 19-9 


DM/70K 


PSA 


CEA 



Prostate Cancer II 
PAP 

PSA 
CR-BB 
Acid Phosphatase 



Breast Cancer 


Testicular Cancer 


Osteoporosis 


LASA-P® 


LASA-P® 


Calcium 




AFP 


Phosphate 


HER2/neu 
in Plasma 


HCG-Beta 


Estradiol 




CA 15-3® 


Progesterone 






ALP 






ALP Isoenzyme 1 






ALP Isoenzyme 2 



A key to the abbreviations used above is provided below: 



AFP: 


Alpha-Fetoprotein 


CA125: 


Cancer Antigen 125 


CA 15-3®** 


Breast Antigens 115D8/DF3 


CA 19-9: 


Carbohydrate Antigen 19-9 


CEA: 


Carcinoembryonic Antigen 


CK-BB: 


Creatinine kinase, BB subtraction 


DM/70K: 


6varian marker NB/70K 


HCG-Beta: 


Human Chorionic Gonadotropin, Beta Sub- 
Unit 


HER 2/neu in 
Plasma: 


c-erb B-2 (HER2/neu) oncoprotein in plasma 


LASA-P®*: 


Lipid-Associated Sialic Acid in Plasma 


PAP: 


Prostatic Acid Phosphatase 


PSA: 


Prostate Specific Antigen 



*LASA-P is a registered trademark of DIANON Systems, Inc. 
**CA 15-3 is a registered trademark of Centocor, Inc. 
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A wide number of diseases may be diagnosed in accordance with the 
method of the present invention. To be suitable for diagnosis by the present 
method, biomarkers for the disease must be quantifiable. The biomarkers 
must also be predictive of the disease and must be statistically significant 
5 relative to one another. The method of the present invention is equally 

suited to the diagnosis of any disease in which serum biomarkers can be 
identified, including but not limited to infectious diseases, and genetic 
abnormalities. 

After determining the serum biomarkers for a disease, the biomarker 

10 values are analyzed by a trained neural network to yield a single diagnostic 

value. The most common neural network architecture for pattern 
classification problems is the feedforward network, which typically consists 
of an input layer, one or more hidden layers, and an output layer. Figs. 1 
and 2 illustrate the arrangement of neurons in two different feedforward 

15 networks. 

The elements that make up each layer of a neural network are 
referred to as neurons or nodes. Inputs are fed forward from the input layer 
to the hidden layers and then to the output layer. The number of neurons in 
each layer is determined before the network is trained. Typically, there is 

20 one input neuron or node for each input variable, and one output node for 

each output The inputs to the neural network are predictor variables. 
These predictor variables can be quantitative or qualitative. Neural networks 
make no data distribution assumptions and can simultaneously use both 
quantitative and qualitative inputs. In the present invention, the biomarker 

25 values, after a linear prescaling to values between 0.0 and 0. 1 , constitute the 

input variables. 

The outputs of the network represent output categories. For 
example, a malignancy may be represented by maximal output of the 
malignant output neuron and silence of the benign neuron, whereas a benign 

30 process is represented by maximal output of the benign neuron and silence 

of the malignant neuron. A simple arithmetic function combines outputs of 
the two neurons to yield a single diagnostic index. In the alternative, a 
single output neuron may be used. An output of greater than 0.5 would 
indicate a malignancy and an output of less than 0.5 would indicate a benign 

35 condition. In this way a diagnostic index is directly obtained. 
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The number of hidden layers and the number of nodes in the hidden 
layers are configurable parameters that have a significant influence on the 
performance of the network. In practice, the optimal number of hidden 
neurons is determined empirically. The means for determining the optimum 
numbers of hidden neurons is well known to those skilled in the art and 
depends on the complexity of the problem being solved. 

In the present invention, the main network model is a multi-layer 
feedforward perceptron using a backpropogation training algorithm. The 
number of hidden layers and the number of neurons in each hidden layer 
was determined to adequately match the level of complexity of the diagnostic 
problem. With the assumption that the samples in the training set are 
representative of all possible situations encountered in real applicationswith 
no significant contradictions, and and the number and stratification of 
samples in the generalization test are statistically adequeate, the criteria 
outlined below is used to determine if a chosen network configuration is 
appropriate. 

If the network continues to fail to correctly classify large portions of 
the samples in the training set, even after many adjustments of training 
algorithm parameters, the network complexity should be increased. 

On the other hand, if the the network acheives a high rate of 
correctly classifying the training set but fails to accurately classify a large 
number of samples in the testing set, network structure is probably too 
complex for the problem being solved, i.e. it has sufficient inherent 
flexibility to fit the training data set, but not sufficient predictive power to 
classify the test data set. If this is the case, the number of neurons in the 
hidden layers should gradually be reduced, or, if there are mutiple hidden 
layers, the hidden layers should be gradually reduced. 

Note that it is not always necessary to have a large training sample 
set. If the samples in a training set have already represented all possible 
cases with adequate statistical significance, the addition of new samples 
genewrally does not increase the amount of information in the training 
samples. Instead it may decrease the useful information to nosie ratio in the 
samples. At the other extreme, too small a training set will generally not be 
able to cover all possible variations in the population. The resultant network 
often simply memorizes all the cases in the training set and does not 
generalize at all. 
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The input and output layers are not directly connected. Every input 
neuron is connected to every neuron in the following hidden layer and 
neuron in ahidden layer is connected to every neuron in the following 
adjacent hidden layer or output layer, depending on the number of hidden 

5 layers. Each of the multiple connections to a particular neuron is weighted. 

In the hidden and output layers, each node sums the input activations, 
multiplied by the respective connection weights plus a bias term. The 
weighted sum then passes through a non-linear output function, typically a 
sigmoidal function, which gives the network the ability to represent complex 

10 non-linear relationships. A neuron fires if the sum of the weighted inputs to 

it are greater than a threshold value. As illustrated in Fig. 3, once a neuron 
is above a threshold, the magnitude of its output is a sigmoid function of the 
net input The end result of activity in the neural network is the net output, a 
complex nonlinear function of the inputs. 

15 In summary and in accordance with the present invention, first the 

values of the biomarkers for a specific disease are determined and scaled. 
The biomarkers are fed forward from the input layer to the hidden layer (or 
layers) and then to the output layer of the neural network. The number of 
neurons in the input layer is determined before the network is trained and 

20 correspond to the number of biomarkers predictive for a specific disease. 

The biomarkers are preselected. There is one input neuron for each 
diagnostic variable or biomarker, and one output neuron for each desired 
output. Other than the identified biomarkers, diagnostic variables may 
include demographic information such as the age, race or sex of the patient 

25 and whether the patient is pre- or postmenopausal. The number of neurons 

in the output layer depends on the type of output desired. The number of 
neurons in the hidden layer is determined euperically during training. 

The neural network used for diagnosing a specific disease must be 
trained to do so. In accordance with the present invention, the neural 

30 network is trained by back propagation. Back propagation refers to the 

technique of training a neural network to accurately model a set of input and 
output examples by determining the best connection weights between the 
values, and is well known in the art. Other techniques which may be used 
to train a neural network for purposes of this invention may include any 

35 other non-linear global optimization technique, such as the genetic search 
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algorithm; however, the feed forward, back propogation network is most 
popular. 

At the beginning of training, the connection weights in the network 
are randomly initialized. The training data is then presented to the network 
one at a time. In accordance with the present invention, the training data 
consists of the biomarker values for a group of patients, and the diagnosis 
for each of those patients. The biomarker values are the input variables used 
to train the network. For each patient, the network uses the patient's 
biomarker values to estimate a diagnosis, which is then compared to the 
actual diagnosis. If the network's diagnosis is correct, then the connection 
strengths and thresholds within the network are not changed, and the next 
patient is presented to the network. If the estimate of the diagnosis is not 
correct, the connection weights and thresholds in both the hidden layer and 
the output layer are adjusted to reduce the size of the classification error. 
After adjustments are made, the next patient is presented. Training proceeds 
until all patients in the training group are correctly classified. If the network 
is not able to completely classify all the data it will train indefinitely, in 
which case training is terminated. 

When training the neural network, the trainer may set the 
decision limits regarding the definition of a classification error, i.e. an 
incorrect diagnosis. The relevant parameter is the error tolerance, which 
specifies how close the estimated output has to be to the actual output to be 
correct. For example, if two output neurons are used and the training 
tolerance is set at 5%, the estimate of malignancy is considered correct if the 
malignant output neuron fires at 95% of maximum and the benign neuron 
fires at 5% of maximum. Similarly, a correct estimate of a benign diagnosis 
means that the benign output neuron firs at 95% of maximum, while the 
malignant neuron fires at 5% of maximum. The methods for determining a 
classification error are well known to those skilled in the art. 

In the preferred embodiment of this invention, if a single output 
neuron is used, a benign normal diagnosis is set at an output of 0.1 and a 
malignant or abnormal diagnosis is set at an output of 0.9. Error tolerance 
is an adjustable parameter and is significant in determining the success of the 
network at making an accurate diagnosis. It is preferred to use either one or 
two output neurons. 
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After the neural network is trained for the desired disease, biomarker 
values from patients with unknown disease conditions are introduced to the 
trained neural network. The neural network then processes the information 
to produce a value corresponding to a diagnosis of the presence or absence 

5 of the particular disease. In accordance with the present invention, this is 

accomplished by using either one single output neuron or two output 
neurons. If two output neurons are used, the output from the two neurons 
are combined to generate a single diagnostic index. 

As illustarted by Fig. 4, in a second embodiment of the present 

10 invention, the diagnostic value obtained by analysis of the biomarkers by a 

trained neural network is further analyzed by a set of heuristic rules in 
combination with additional patient information. The additional patient 
information includes things such as family medical history and demographic 
information. This data is then processed to yield a second single diagnostic 

15 value. 

This invention is further illustrated by the following examples, 
which are not to be construed in any way as imposing limitations upon the 
scope thereof. On the contrary, it is to be clearly understood that resort ma> 
be had to various other embodiments, modifications, and equivalents 

20 thereof, which, after reading the description herein, may suggest themselves 

to those skilled in the art without departing from the spirit of the present 
invention and/or the scope of the appended claims. 

In the following examples which utilize a neural network in the 
analysis of the data, a Neuralshell 2, Release 1.5 (Ward Systems Group, 

25 Inc.) neural network development program was used for the training of the 

neural network on a Pentium 60 mhz computer (Magitronic, Inc.). 

1 



The following example describes the training of a neural network to 
30 diagnose colon cancer. 

A total of 167 patients ( 48 with diagnosed colon cancer, and 119 
normal) were divided into 2 groups, a training set and a generalization 
testing set. The training set contained approximately 60% of the patients (90 

♦ 

patients, 27 with diagnosed colon cance, and 63 normal), and the 
35 generalization testing set contained approximately 40% of the patients (77 

patients, 23 with diagnosed colon cancer, and 54 normal). 
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The initial network architecture was selected based on the level of 
complexity of the classification task. A multi-layer feedforward network 
was used. Selection of the initial architecture involved the selection of the 
number of hidden layers and the number of neurons in each hidden layer. 
Several trial iterations were performed to determine and adequeate 
configuration that showed good results on both the training sample set and 
the generalization test sample set. This particular neural network had 5 input 
neurons, one hidden layer having 10 neurons and one output neuron. 

Initially, connection weights among the neurons were randomly set 
The neural network had five input neurons, corresponding to five input 
variables for diagnosing colon cancer: LASA-P®; CA 19-9, CEA, Age, 
and Sex. During training, the five input variables for each patient were first 
linearly scaled into the continuos range between 0.0 and 1 .0. The resultant 
five numbers were then presented as an input vector to the input neurons of 
the artificial neural network. 

For each of the input vectors, the network generated an output based 
on the connection weights among the network neurons. The output can be a 
single value or a vector of numbers, depending on the number of output 
neurons used. Each neuron in the network participated in the output 
calculation by passing the sum of all inputs to the neuron through a non- 
linear s-shaped function (often a ligistic function) and sending the result to 
each and everyone of the neurons in the following adjacent layer. The 
generated output was compared to the desired "target" output. A value of 
0.1 corresponded to a diagnosis of normal and an output of 0.9 
corresponded to a diagnosis of of abnormal. The difference was used to 
calculate an error term to guide the training algorithm, i.e., the 
backpropogation algorithm, in the adjustment of network connection 
weights in an attempt to reduce the differences between network outputs and 
target values over the training sample set. 

To counteract over learning, which can occur when a neural network 
is too large for the problem and trains too long, the network was tested on 
the generalization test set every 200 iterations. The networks for which 
generalization testing showed reduced root mean squared error over 
previous tests were saved and used in later calculations or generalizations. 
This procedure was repeated until generalization failed to improve over 2000 
successive tests. 
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10 



To accelerate the rate of solution during the initial stage of training, 
the learning rate for the hidden layer was initially set at 0. 1 and momentum 
was set at 0.05, for the first several iterations. After the first several 
iterations, the learning rate was lowered to 0.01 and momentum to 0.0005 
to narrow in on the appropriate solution. 

After training, the neural network correctly identified 95% of the 
patients without cancer in the training data set, and 92% of the patients with 
cancer. 

When presented with the generalization test results, the trained 
neural network correctly identified 100% of the people with normal 
diagnosis and 95% of the people with cancer diagnosis. 



15 



20 



25 



Example I 

The following example illustrates the value of the present invention 
for diagnosing ovarian cancer. 

The biomarker values for LASA-P®, CA 125, and DM/70K were 
measured in a patient suspected of having ovarian cancer. Two of the 
biomarker values fell within the normal reference range, and one value fell 
outside the normal reference range. The values were scaled and analyzed by 
a trained neural network, yielding a diagnostic index of 75, indicating the 
presence of ovarian cancer. The diagnotic index generated by the neural 
network (between 0.1 and 0.9) is multiplied by 100 to obtain a whole 
number. The normal reference range for the diagnostic index was 0-50. 
The biomarker values measured, together with the normal reference range 
for those values are provided below: 



Biomarker 


Patient Value 


Normal Reference 






Range 


LASA-P® 


18 


15-20 mg/dl 


CA 125 


16 


0-35 U/ml 


DM/70K 


45 


0-35 U/ml 


Neural Net 
Diagnostic Index 


75 "j 


0-50 



30 



The patient underwent surgery to treat, the cancer. Fig. 5 graphically 
illustrates the change in the neural network diagnostic index over time, 
before and after surgery. 
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Example 3 

The following example illustrates the value of the present invention 

for diagnosing prostate cancer. 

The biomarker values for LASA-P®, PAP, and PSA were measured 
in a patient suspected of having prostate cancer. All three of the biomarker 
values fell within the normal reference range. The values were scaled and 
analyzed by a trained neural network, yielding a diagnostic index of 90, 
indicating the presence of prostate cancer. The normal reference range for 
the diagnostic index was 0-50. The biomarker values measured, together 
with the normal reference range for those values are provided below 




The patient underwent surgery to treat the cancer. Fig. 6 graphically 
illustrates the change in the neural network diagnostic index over time, 
before and after surgery. 

Enamels. 4 

The following example illustrates the value of the present invention 

for diagnosing colon cancer. 

The biomarker values for LASA-P®, CA 19-9, and CEA were 
measured in a patient suspected of having colon cancer. All three of the 
biomarker values fell within the normal reference range. The values were 
scaled and analyzed by a trained neural network, yielding a diagnostic index 
of 93, indicating the presence of colon cancer. The normal reference range 
for the diagnostic index was 0-50. The biomarker values measured, 
together with the normal reference range for those values are provided 
below: 
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Biomarker 


Patient Value 


Normal Reference Range 


LASA-P® 


17.9 


15-20 mg/dl 


CA19-9 


5.0 


0-37 U/ml 


CEA 


2.3 


0-2.5 ng/ml 


[ Neural Net 

| Diagnostic Index 




| 0-50 



The patient underwent surgery to treat the cancer. Fig. 7 graphically 
illustrates the change in the neural network diagnostic index over time, 
before and after surgery. 

5 

Example g 

The following example illustrates the value of the present invention 
for diagnosing breast cancer. 

The biomarker values for LASA-P®, CEA, HER2/neu in Plasma 

10 and CA 15-3® were measured in a patient suspected of having breast 

cancer. All four of the biomarker values fell within the normal reference 
range, although two of the biomarker values were on the high end of the 
normal range. The values were scaled and analyzed by a trained neural 
network, yielding a diagnostic index of 91, indicating the presence of breast 

15 cancer. The normal reference range for the diagnostic index was 0-50. The 

biomarker values measured, together with the normal reference range for 
those values are provided below: 



Biomarker 


Patient Value 


Normal Reference Range 


LASA-P® 


18 


15-20 mg/dl 


CEA 


1.0 


0-2.5 ng/ml 


HER2/neu in Plasma 


19 


0-20 U/ml 


CA 15-3® 


8 


0-30 U/ml 


Neural Net 

Diagnostic Index | 


1 91 


0-50 



20 The patient underwent surgery to treat the cancer. Fig. 8 graphically 

illustrates the change in the neural network diagnostic index overtime, 
before and after surgery. 
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Example 6 

This example illustrates the construction and training of a neural 
network for diagnosis of ovarian cancer. 

Table I below provides the data used to train the network (age, 
LASA-P®, CA 125, and DM/70K), the actual histological diagnosis, and 
diagnostic index determined by the neural network. A histological index of 
0.1 corresponds to a normal diagnosis and a histological index of 0.9 
corresponds to an abnormal diagnosis. A neural net diagnostic index less 
than .5 corresponds to a normal diagnosis. A neural net diagnostic index 
greater than .5 corresponds to an abnormal diagnosis. 

TABLE I 

Training Data for Ovarian Cancer Diagnostic System 
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84 

S3 

89 

92 

2042 
2043 
2044 
2046 
2047 
2048 
2049 



33 

26 

26 

25 

68 

63 

84 

72 

61 

44 

60 



17.4 
15.6 
20.7 
20.7 
16.4 
14.5 
17.6 
17.8 
13.4 
20.8 
17.8 



19 



21 



12 



22 



23 

8 

146 

21 

29 



0_ 

0 



0_ 

17 

15 

15 

0 

0 



0.1 

0.1 

0.1 

0.1 

0.9 

0.9 

0.9 

0.9 

0.9 

0.9 

0.9 



0.305 
0.156 
0.205 
0.164 
0.852 
0.733 
0.931 
0.813 
0.852 
0.567 
0.822 



An evaluation of the training data set (n=61) yielded a sensitivity of 96<7, 
and a specificity of 100%. 

Table II provides the data used to test the network trained with the 

data provided in Table I. 



TABLE II 
Testing Data for Ovarian Cancer 



ic System 



Patient 
Number 


Age 


^ASA-Pd 


CA 125 1 


DM/70K 


Histological 
Diagnosis 1 


Neural Net 1 
Diagnostic 1 
Index 1 
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43 


14.1 1 


8 1 


J_5 


0.1 


0.371 1 


1104 


28 


15.9 


~7l 


0 


OA 


0.211 1 


1 105 


46 


15.9 1 


15 1 


0 
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1 0.529 1 


1 107 


44 


23.6 1 


10 . 1 
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0.1 I 


1 0.416 1 
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1 9 
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1° 
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0.9 


|| 0.823 
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The testing data (n=59) yielded a sensitivity of 92% and specificity 
of 91%. 

The combined results of the training and testing data yielded a 
5 sensitivity of 94% and specificity of 96%. 

Fig. 20 provides the equations describing the weight matrices 
between the layers in the neural network 

Table III below compares the sensitivity and specificity of a 
diagnosis obtained by the method of the present invention with the 
LO sensitivity and specificity of each individual biomarker. As is evident from 

the results, the individual biomarkers may be combined by this method to 
yield a single diagnostic index having greater sensitivity than the individual 
tests for purposes of diagnosis. 
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15 



20 
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TABLE in 

Comparision of Ovarian Cancer Diagnostic Methods 



BioMarker Test 


Sensitivity 


Specificity 


Carbohydrate Antigen 125 
(CA 125) 


74%* 


99%* 


Dianon Marker 70K 
(DM/70K) 


68%* 


99%* 


Lipid Associated Sialic Acid 
(LASA-P®) 


79%* 


96%* 


Neural Network Diagnostic 
Index for Ovarian Cancer 


94%** 


96%** 



* BioStatislics provided by Dianon Systems, Inc. 
** BioStatistics based on initial pilot study performed by Horus 
Therapeutics, Inc. 

Example 7 

This example illustrates the construction and training of a neural 
network for diagnosis of osteoporosis. 

Fig. 9 provides the data used to develop, i.e. train and test, the 
neural network to diagnose osteoporosis. The biomarkers selected included 
age, calcium, phosphate, estradial (ETWO), progesterone, total alkaline 
phosphatase, total intestinal alkaline phosphatase, and % liver alkaline 
phosphatase. Fig. 9 further includes the diagnostic index obtained by the 
neural network. 

Fig. 10 provides the data use to test the network trained with the data 
in Fig. 10, and the neural network diagnostic index obtained. 

The second embodiment of the present invention also relates to a 
method for diagnosing diseases in humans or animals using blood 
concentrations of specific, predetermined blood constituents. More 
specifically, the present invention relates to a method for detecting the 
occurrence of osteopenia, diagnosing osteoporosis and determining its 
severity and underlying cause. The present invention also facilitates the 
periodic monitoring of specific physiological functions which may indicate 
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the onset of osteopenia and correlates to bone mineral density measurements 
determined by various standard methods. 

In practicing one aspect of the present invention, the severity of 
disease in a set of humans or animals with varying severity of disease is 
5 measured by a standard method or methods. The measurement is then 

assigned a numerical value corresponding to a severity scale. The scale 
ranges from humans or animals with no disease, to humans or animals with 
severe disease. The scale is preferably a numerical scale. For example, one 
could assign a value which corresponds to normal or slight disease, another 
10 value which corresponds to moderate disease and a third value which 

corresponds to severe disease. 

The concentration of a predetermined set of blood constituents in the 
set of humans or animals with varying severity of disease is then 
determined. According to the present invention, it is preferable to measure 
15 the blood constituents in the same set of humans or animals in which the 

severity of disease was measured by the conventional method or methods. 

After determining the blood concentrations of the predetermined set 
of blood constituents in the set of humans or animals with varying severity 
of the disease, mathematical manipulations are performed in which the 
20 numerical value obtained using the conventional diagnosis is set to equal the 

result of a mathematical model utilizing the set of blood constituent 
measurements. For example, the relationship can be made using a multiple 
linear regression analysis. It is to be understood that other mathematical 
models could be used to determine a correlation between the serum 
25 constituent concentrations and a standard method of determining severity of 

the disease. The concept of using mathematical models to determine the 
relationship or correlation is considered to be part of the present invention. 
Standard statistical analyses that are well-known to those of ordinary skill in 
the art are then performed to determined the confidence level of the 
30 correlation between the diagnosis by conventional means and the set of 

blood constituents. These statistical analyses include chi-square tests. 

An example of practicing one embodiment of the present invention is 
a method for diagnosing osteopenia in a human or animal. The method 
preferably utilizes six blood constituents. These constituents are calcium, 
35 phosphate, total alkaline phosphatase, an alkaline phosphatase isoenzyme, 

estradiol, and progesterone. The alkaline phosphatase isoenzymes preferred 
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for practicing the present invention include lymphocyte-derived alkaline 
phosphatase isoenzyme and bone, liver or intestinal alkaline phosphatase 
isoenzymes. The present invention includes calculating a bone density 
quotient using the aforementioned six blood constituents by entering the 

5 values for the tests into an algorithm that is calculated using a multiple linear 

regression analysis. Optionally, the age of the patient may also be 
incorporated into the equation. The bone density coefficient derived from 
the algorithm allows one to diagnose the osteopenic state of the patient, 
including the severity of the disease. 

10 In addition to diagnosing the osteopenic state of the human or 

animal, an indication of the underlying cause of the osteopenia can be 
determined using the present invention. For example, by practicing the 
present invention as described herein, one can determine whether the 
osteopenia in a human or animal is caused by post-menopausal lack of 

15 estrogen or is caused by some other condition, such as cancer. This allows 

the attending physician to be better able to prescribe the appropriate 
treatment for the osteopenia. 

Five of the serum tests that are used in the present invention are tests 
that are commonly performed by clinical laboratories. The test for t- 

20 lymphocyte derived alkaline phosphatase is experimental only; however, the 

test for blood, liver and intestinal alkaline phosphatase isoenzymes are also 
known. The type of test used to determine the six serum constituents is not 
critical to the present invention as long as the tests give accurate blood 
concentrations of the constituents being measured. 

25 

Calcium & Phosphorus: 

Maintenance of calcium and phosphorus homeostasis involves the 
participation of 3 major organs: the small intestine, the kidney, and the 
skeleton, and is regulated by various hormones. Calcium enters the body 

30 through the diet and is absorbed into the circulation from the small intestine. 

Calcium absorption occurs by two processes, active transport and passive 
transport. Approximately 98% of the calcium and 85% of the phosphorus 
in the adult is present in the skeleton primarily as hydroxyapatite, which is a 
crystal lattice compound of calcium, phosphorus, and hydroxide. 

35 Remaining calcium is present in extra cellular fluid, some types of tissue, 

and skeletal muscle. Phosphorus is combined with lipids, proteins, 
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carbohydrates, and other organic substances. Of critical importance to 
calcium homeostasis is the fact that less than 1% of total skeletal reservoir of 
calcium is rapidly exchangeable with extracellular fluid. In addition to its 
obvious importance in skeletal mineralization, it is also vital for blood 
5 coagulation, neuromolecular condition, maintenance of normal tone and 

excitability of skeletal and cardiac muscle, and preservation of all membrane 
integrity and permeability, particularly in terms of sodium and potassium 
exchange. 

Free or ionized calcium accounts for 50% of total calcium. About 
10 5% of total calcium is complexed with a variety of anions, particularly 

phosphate and citrate. The remaining 45% of calcium is bound to plasma 
proteins. Both ionized calcium and calcium complexes are freely dialyzable. 
Acidotic and alkalotic conditions adversely affect the ionized calcium level in 
the blood. In metabolic bone disease such as hyper- or 
15 hypoparathyroidism, Pagets of bone, Vitamin D deficiency, and renal 

osteodystrophy, calcium and/or phosphorus levels are significantly altered. 

Serum Calcium Determination: 

In one method of practicing the present invention, the procedure for 

20 the determination of calcium is based on the interaction of the chromogenic 

agent o-cresolphthalein complexone (Sigma Diagnostics Calcium Agent, 
Sigma Chemical Co., St Louis, MO ) which complexes with calcium cation 
in an alkaline medium to form a purple colored complex which has an 
absorbance maximum at 575 nm. The intensity of the color measured at 575 

2 5 nm is directly proportional to the calcium concentration in the given sample. 

The o-cresolphthalein complexone contains 8-hydroxyquinoline which 
prevents interference from magnesium ions. Robertson, W.G., et al., 
"Calcium Measurements In Serum and Plasma - Total and Inonized," Crit. 
Rev. Clin. Lab. Sci., Vol. 11, pg. 271, (1979); Sigma Diagnostics Product 

30 Insert: Calcium, Procedure No. 587, Reissued May 1989. The calcium 

level in serum is in the range of about 9.2 to 1 1.0 mg/dl which is composed 
of 3 distinct fractions. 
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Specimen Collection: 

In one method of practicing the present invention, serum or 
heparinized plasma are suitable for the test. Anticoagulants other than 
heparin should not be used. 

5 

Procedure: 

Sigma Diagnostics Calcium Reagent is used with Abbott Spectrum 
High Performance Diagnostic System for the quantitation of calcium in 
serum. It is generally considered good laboratory practice to run a 
10 calibration, linearity, and a quality control prior to the determination of 

patient samples. Calcium levels up to 16 mgs % can be measured by this 
method. 

Serum Phosphorus Determination: 

15 In one method of practicing the present invention, the procedure for 

the determination of serum phosphorus is based on the interaction of 
inorganic phosphorus with ammonium molybdate in the presence of sulfuric 
acid. Barnhill S., etal, "Osteoporosis: A Possible Autoimmune Etiology," 
Ann. of Clin. Lab. ScL, Vol. 17, pg. 255, (1987). This reaction produces 

20 unreduced phosphomolybdate complex. The absorbance of this complex at 

340 nm is directly proportional to the inorganic phosphorus present in the 
given sample. (Sigma Diagnostics Product Insert: Phosphorus, Inorganic, 
Procedure No. 360-UV. Previous Revision December 1985. Reissued 
July (1988)). Most of the phosphorus in extracellular fluid is inorganic and 

25 in the range of 2.4 to 4.7 mg/dL 



Specimen Collection: 

Serum or heparinized plasma are preferred for the test 
Anticoagulants other than heparin should not be used. 
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Procedure: 

Sigma Chemical Company's phosphorus inorganic reagent can be 
used with the Abbott Spectrum High Performance Diagnostic System for 
qualification of the inorganic phosphorus in serum. It is customary to 
perform calibration, linearity, and quality control studies prior to the 
determination of test samples. Phosphorus levels up to 12 mgs % can be 
measured by this procedure. 

Total Alkaline Phosphatase: 

The demonstration that bone is rich in alkaline phosphatase (ALP) 
and that normal serum contains the same or a similar enzyme led to the study 
of serum ALP levels in patients with disease of bone. See, Courpron P. , 
"Bone Tissue Mechanisms Underlying Osteoporosis," Orthop. Clin. North 
Am., Vol. 12, page 513, (1981). This is especially true in osteitis 
deformans, hyperparathyroidism and bone neoplasm. ALP is also elevated 
in hepatic disease; however, this can be distinguished by other corroborative 
laboratory procedures and clinical features. In some situations, as in 
osteoporosis, ALP per se may not be above the reference range; however, 
the isoenzyme of ALP is increased. It is well understood that total serum 
ALP in normal subjects consists of isoenzymes contributed from liver, 
bone, renal, pulmonary, placental and intestinal sources, among others. 
The isoenzymes and isoforms of various tissue origins can be further 
separated and visualized by an Isoelectric Focusing Technique. Estimation 
of total serum ALP by isoelectric focusing techniques offers great potential 
in the investigation of metabolic bone diseases. 

Serum Alkaline Phosphatase Determination: 

In one method of practicing the present invention, serum ALP 
activity can be measured using various phosphate esters as substrates. 
Sigma Chemical Company's alkaline phosphatase reagent measures serum 
ALP activity by a kinetic method. The reagent for the test contains p- 
nitrophenyl phosphate, carbonate buffer, magnesium ions and mannitol. 
Mannitol present in the reagent acts as a phosphate acceptor during the 
enzyme reaction. McComb R.B., et al., Alkaline Phosphatase. Plenum 

Press, New York, (1979); Gundbery, Af., AlKnlinr Phosphatase and 

Osteocalcin. Primer on M etahnlie Bop» DjgfiaSSS anri Disorders of Mineral 
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Metabolism 1st Edition, Published by Am. Soc. For Bone and Mineral Res. 
Editor-Murry J. Favus, pgs. 74-76, (1990); Sigma Diagnostics Product 
Insert: Alkaline Phosphatase (ALP), Procedure No. 245, Reissued March 
1989; Kaplan TA.,et al, Clinical Chemistry . St. Louis, e.v. Mosby 
5 Company, (1987). 

Specimen Collection: 

Preferably, serum or heparinized plasma is used for the total alkaline 
phosphatase and alkaline phosphatase isoenzyme determination tests. 
10 EDTA, oxalate, citrate and fluoride are inhibitors of ALP and are not 

suitable anticoagulants. 



Procedure: 

Alkaline phosphatase (ALP) diagnostic reagent manufactured by 
15 Sigma Chemical Company is used with the Abbott Spectrum High 

Performance Diagnostic System for the quantitation of inorganic 
phosphorus in the serum. Serum ALP hydrolyzes p-nitrophenyl phosphate 
to p-nitrophenol and inorganic phosphate. The hydrolysis occurs at alkaline 
pH and the absorbance is directly proportional to the ALP activity of the 
20 serum sample. As is the usual practice, calibration, linearity and quality 

control assays should be performed prior to the determination of test 
samples. 

ALP levels up to 1200 U/L can be measured by this procedure. 
Normal ranges are: 
25 Infants 50-165 (U/L) 

Adult 20-70 (U/L) 

Child 20-150 (U/L) 

60 years 30-75 (U/L) 



30 Alkaline Phosphatase Isoenzyme: 

Alkaline phosphatase isoenzyme levels have been primarily used to 
aid in the differential diagnosis of liver and bone disorders and also to 
indicate placental growth patterns. Isoelectric focusing of alkaline 
phosphatase isoenzyme on a precast gel gives 12 discernible bands of 
35 different organ and tissue origin. Of these 12 Bands, Band 10, appearing at 

4.73 isoelectric point, is of placental origin. It is released by placental 
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syncytiotrophoblast and/or macrophage, and is T-lymphocyte mediated. In 
the alternative, samples for ALP isoenzyme can be analyzed by 
electrophoresis methods which identify mainly liver, bone, and intestinal 
isoenzymes. The Ciba Corning Alkaline Phosphatase Isoenzyme (ALP) 
5 System may be used for the qualitative and quantitative determination of 

ALP isoenzymes in human serum by electrophoresis. 

The liver isoenzyme is the most frequently encountered alkaline 
phosphatase isoenzyme in pathologically increased serum alkaline 
phosphatase activity. Increased liver isoenzyme is encountered in a variety 
10 of liver and hepatobiliary diseases. These include hepatic cirrhosis, primary 

biliary cirrhosis, congestive cirrhosis, and hepatic carcinoma. Elevation of 
the liver isoenzyme is also a sensitive indicator of cholestasis and liver 
infiltration. 

Another of the alkaline phosphatase isoenzymes that can be used in 
1 5 practicing the present invention is lymphocyte-derived alkaline phosphatase 

The preferred method of measuring this isoenzyme is by isoelectric focusing 
electrophoresis. A commercially available isoelectric focusing apparatus that 
is capable of measuring alkaline phosphatase isoenzymes is made by 
IsoLab, Inc. (Resolve® -ALP, IsoLab, Inc., Akron, Ohio). On the 
20 Resolve®-ALP isoelectric focusing apparatus, the alkaline phosphatase 

isoenzyme that is used in the present invention resolves as band 10 in the 
electrophoretic pattern. In the preferred bone density algorithm defined 
hereinbelow, the value for the band 10 alkaline phosphatase isoenzyme was 
assigned 0 if the band were missing or very weak, and 1 if the band were 
25 present. The electrophoresis gels can also be scanned with a densitometer 

and more quantitative values can be assigned to the alkaline phosphatase 
isoenzyme concentration. As used herein, the band 10 alkaline phosphatase 
isoenzyme shall be designated lymphocyte-derived alkaline phosphatase 

isoenzyme, or I-Alkp. 

It is to be understood that other methods of measuring I-Alkp could 
be used. These methods include, but are not limited to, enzyme-linked 
immunoassay techniques (ELISA), radioimmunoassay techniques, affinity 
columns, and isoelectric focusing columns. It is important to note that 
measurement of band 10 on the Resolve®-ALP apparatus does not correlate 
at all with osteopenia in general. There are many other abnormal conditions 
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which result in a higher than normal band 10 lymphocyte-derived alkaline 
phosphatase isoenzyme (see Barnhill, etaL and Griffith, et ai, supra.). 

A number of pathological conditions can lead to an elevated level of 
bone alkaline phosphatase isoenzyme. The highest levels of bone 

5 isoenzymes are usually found in Pagets disease. Increased levels are also 

encountered in rickets, bone cancer, osteomacacia and celiac sprue. Renal 
disorders can also result in increased levels. These include renal failure, 
primary hyperparathyroidism, secondary hyperparathyroidism induced by 
long term hemodialysis and malabsorption. 

10 Increased levels of intestinal alkaline phosphatase are encountered in 

a variety of diseases of the digestive tract These include intestinal infection 
and ulcerative lesions of the stomach, duodenum, small intestine and colon. 



Procedure: 

15 In the preferred method of practicing the present invention, the 

alkaline phosphatase isoenzymes are separated by electrophoresis in a 
buffered agarose system. After electrophoresis the isoenzymes are detected 
by incubating the gel with a fluorescent compound such as 4- 
methylumbelliferyl phosphate. The fluorescence formed during the 

20 interaction is quantitated using a Ciba-corning 710 densitometer at 385 nm 

Total alkaline phosphatase activity of the patient is required to interpret the 
liver, bone and intestinal ALP isoenzymes. 



Estrogen and Progesterone: 

25 The estrogens are steroids that have a ring containing three 

unsaturated double bonds. The ovary, as well as the testes and adrenal 
gland, has the capacity to synthesize estrogens from androgens, 
androstenedione and testosterone. During the follicular phase of the 
menstrual cycle, ovarian secretion represents only one third of total estrogen 

30 production. In contrast to estradiol, which is secreted almost entirely by the 

ovary, most estrone is derived from peripheral conversion of 
androstenedione and from estradiol metabolism. During menopause, 
estradiol concentrations steadily decrease to approximately 15 percent of 
premenopausal levels. In healthy postmenopausal women, the ovaries do 

35 not secrete significant quantities of estrogens, and virtually all estrogen 
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produced is from peripheral conversion of androstenedione made by the 
adrenal. Estradiol is present in the serum as follows: 



Follicular Phase 




Early 


30-100 ng/L 


Late 


100-400 ng/L 


Luteal Phase 


50-150 ng/L 


Postmenapausal 


<20 ng/L 



10 



15 



In menstruating females, progesterone is secreted mainly by the 
corpus luteum of the ovary. It is partially responsible for cyclic changes in 
the endometrium that are necessary for attachment and growth of an 
embryo. Progesterone levels are low prior to the mid-cycle gonadotropin 
surge. Shortly after the gonadotropin surge, they begin to rise rapidly, 
reaching peak levels during the middle of the luteal phase. Thereafter, a 
progressive fall occurs with barely detectable progesterone levels prior to 
menses. Although progesterone in large amounts produces a negative 
feedback on gonadotropin secretion, it is not the major component in the 
negative feedback system of ovarian steroids. Function of the corpus 
luteum can be assessed by measuring serum progesterone concentration. 
Progesterone is present in the serum as follows: 



Follicular Phase 


0.1 - 1.5 ng/L 


Luteal Phase 


2.5 - 28.1 ng/L 


Mid-Luteal Phase 


5.7 - 28.1 ng 


Over 60 Years 


0.0 - 0.2 ng/L 



20 



Procedure: 

In a preferred method of practicing the present invention, progesterone 
measurements can be obtained by a radioimmuno-assay method using 
antibody coated tubes (Diagnostic Products, Los Angles, California). 
Estradiol measurements can be performed by a Microparticle Enzyme 
Immunoassay (IMX), available from Abbott Diagnostics, Abbott Park, IL. 
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Bone Mineral Density (BMD) Measurements 

The conventional methods of diagnosing osteopenia which may be 
used when practicing the present invention are outlined below: 

Dual-Photon absorptometry (DP A) is widely used to assess bone 

5 mineral content and bone mineral density. Johnston Conrad C, et al. 9 

"Clinical Use of Bone Densitometry", New Eng. J. of Med, Vol. 324, pgs. 
1 105-1 109, (1991). DPA uses transmission scanning with photons from a 
radioisotope source, such as 153 Gd, that emits two energy peaks, thus 
allowing bone density to be measured independent of soft tissue. DPA 

10 measurements are performed on lumbar spine (L1-L4) and femur (femoral 

neck, Ward's Triangle and trochanter) and the average determined 
separately. The overall average of both hip and spine can also be 
determined. BMD can be measured by DPA using - Gadolinium as the 
source (Lunar DP3, by Lunar Radiation Corporation, Madison, Wisconsin). 

15 The present method employs transmission scanning using 44 and 

100 KeV photon energies from a one Ci 153 Gd source to allow computation 
of the mineral content of bone independent of soft tissue thickness. Bone 
mineral density, expressed in g/cm2, is derived by dividing bone mineral 
content (BMC) by the projected area of the scanned bone. Peppier W.W., 

20 et a/., "Total Body Bone Mineral and Lean Body Mass by Dual-Photon 

Absorptometry: 1. Theory and Measurement Procedure," CALCIP, Issue 
Int., Vol. 33, pg. 353, (1981); Shipp C, et aL, "Precision of Dual-Photon 
Absorptometry," CALCIP, Issue Int., Vol. 42, pgs. 287-292, (1988). 
During spine scans, the detector moves in a rectilinear pattern at a rate of 

25 5mm/sec and with scan lines 4.5 mm apart. The BMC and BMD are 

calculated in lumbar vertebrae 1 through 4 (including intervertebral discs) 
with a software version supplied by Lunar Radiation Corporation of 

w 

Madison, Wisconsin. 

During femur measurements, the scanner moves at a rate of 2.5 

30 mm/sec and a step distance of 2.5 mm. The BMC and BMD of the neck, 

Ward's Triangle and trochanteric regions of the proximal femur are 
calculated using a femur software version supplied by Lunar Radiation 
Corporation. The femoral neck region of interest (ROI) is that band about 
1 .5 cm wide across the neck of the bone perpendicular to the neutral axis 

35 with the lowest density. Ward's Triangle is defined as a square ROI (about 

1.5 x 1.5 cm) with the lowest density within the proximal femur region. 
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Ward's Triangle is predominantly trabecular bone and contains the least 
amount of bone mineral within the neck region. Carter D.R., et al., 
"Relationship Between Loading History and Femoral Cancellous Bone 
Architecture," J. Biomechanics, Vol. 22, pgs. 231-244, (1989). In the 
proximal femur, the ROI is usually the area 1.5 cm wide across the entire 
femoral neck. Additional regions are defined by the software in the lower 
density Ward's Triangle region, and in the region of the greater trochanter. 
Bone loss in the proximal femur begins in the Ward's Triangle region and 
proceeds outward from there. {Brown D., et al, "Mechanical Property 
Distributions in the Cancellous Bone of the Human Proximal Femur," Act. 
Orthop. Scand., Vol. 51, pgs. 429-437, (1990). This makes the region an 
early indicator of bone loss, but the higher variance in measuring it, 
compared to the neck region, makes the latter zone, a better discriminator. 
However, Ward's Triangle is the lowest density area at the point where the 
neck and greater trochanter meet, a primary hip fracture site. This 
operational ROI may not correspond exactly to the anatomic Ward's 
Triangle region but does provide a repeatable measurement. The width of 
the neck ROI and the size of the Ward's Triangle ROI are actually 
proportional to the measured size of the femoral neck. It has been shown 
that the density of the Ward's Triangle area is substantially reduced in hip 
fracture patients compared with age matched controls. (Vose G. P., et al, 
"Femoral Neck Fracturing its Relationship to Radiographic Bone Density," 
J. Gerontol, Vol. 20, pgs. 300-305, (1965). DPA measurements have a 
precision of 1-3% and the scan can be completed in about 20 minutes. Bone 
density of vertebrae correlates well with risk for vertebral fracture, and bone 
density of areas in the proximal femur correlates well with risk for hip 
fracture. DPA has a low radiation dose (<10 mrem to skin, and 2 mrem to 
marrow). DPA has approximate sensitivity of spine and hip of 50% and 
53%, respectively, at 95% specificity (i.e. % of fracture cases below 5th 
percentile). Moreover, DPA is not subject to systematic errors introduced 
by variable osteoid and variable marrow. 

Bone density can also be measured by radiographic absorptometry. 
Radiographic absorptometry is a method of measuring bone density which 
is well known to those of ordinary skill in the art. Other methods of 
measuring bone density include quantitative computed tomography and 
direct measurement of bone density. 
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Direct measurement of bone density can be obtained by measuring 
the bone density in cadavers. Thus, another method of determining the 
correlation between the blood concentrations of certain predetermined blood 
constituents and bone density is to measure the bone density of a set of 
5 cadavers from which blood concentrations of predetermined blood 

constituents are known. The multiple linear regression analysis can then be 
performed and a bone density coefficient relationship can easily be 
determined. 

In correlating the bone density measurements to concentration of 
10 blood constituents, blood concentrations of calcium, phosphate, total 

alkaline phosphatase, an alkaline phosphatase isoenzyme estradiol, and 
progesterone are measured. Liver, bone and intestinal isoenzymes can be 
used. Band 10 alkaline phosphatase isoenzyme. A mathematical 
relationship between the concentrations of blood constituents and bone 
15 density, as measured by radiographic absorptometry or other standard 

method of measuring bone density, is determined by performing a multiple 
regression analysis with the following model: 

Bone density coefficient =fc 0 + MCa + Z?2P + &3E2 + £>4Pg +*5Alkp + 
20 £>6l-Alkp + *>7(CaCa) + &8(PP) + i>9(E2E2) + £>lO(Pg + p g) + 

fell(AlkpAlkp) + *>12(CaP) + &i3(CaE2) + £>14(CaPg) +*>i5(CaAlkp) + 
frl6(CaI-Alkp) + fcl7(PE2) + i>18(PPg) + &19(PAlkp) + Z>20(PI-Alkp) + 
£?2l(E2Pg) + £>22(E2Alkp) + *>23(E2I-Alkp) + £>24(PgAlkp) + &25(PgI- 
Alkp) + fc26(AlkpI-Alkp) 

25 where: 

Ca = Serum concentration of calcium 

P = Serum concentration of phosphorus 

E2 = Serum concentration of estradiol 

PG = Serum concentration of progesterone 
30 Alkp = Serum concentration of total alkaline phosphatase 

I-Alkp = Serum concentration of alkaline phosphatase isoenzyme. 

Calculation of the mathematical model utilized the Systat® statistical 

package (Systat: Inc., Evanston, IL). The multiple linear regression 

analysis is an iterative process which calculates the correct coefficients so 
35 that the result of the algorithm using the blood concentrations correlates to a 
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high degree with the result of the bone density measurement by radiographic 
absorptometry. 

For diagnosing osteopenia, the general form of the preferred 
algorithm that is used in the present invention is as follows: 

Bone density coefficient = b 0 + fciCa + b2P + &3E2 + J?4Pg +*>5Alkp + 
Z?6l-Alkp + &7(CaCa) + i>8(PP) + &9(E2E2) + &10(Pg + Pg) + 
fcll(AlkpAlkp) + fcl2(CaP) + M3(CaE2) + fci4(CaPg) +£>i5(CaAlkp) + 
fcl6(CaI-Alkp) + fcl7(PE2) + i>18(PPg) + *19(PAlkp) + &20(PI-Alkp) + 
f?2l(E2Pg) + ^22(E2Alkp) + £>23(E2I-Alkp) + £>24(PgAlkp) + *?25(PgI- 
Alkp) + i?26(AlkpI-Alkp) 
where: 

bo is a constant 

Ca = Serum concentration of calcium 

P = Serum concentration of phosphorus 

E2 = Serum concentration of estradiol 

PG = Serum concentration of progesterone 

Alkp = Serum concentration of total alkaline phosphatase 

I-Alkp = Serum concentration of alkaline phosphatase isoenzyme. 

It is to be understood that concentration of I- Alkp can be replaced 
with serum concentration of liver, bone or intestinal derived alkaline 

phosphatase isoenzymes. 

In practicing the present invention, the values obtained from the 
patient for the six blood constituent concentrations are inserted into the 
algorithm as indicated. The mathematical manipulation is performed. The 
resulting number is called a bone density coefficient and is then placed in the 
severity scale. This results in a bone density probability quotient which 
correlates to a high degree to bone density as measured by radiographic 
methods. 

As seen in the algorithm that is considered part of the present 
invention, there are a number of coefficients which are a part of the 
algorithm. It is to be understood that these coefficients can change if the 
bone density coefficient is correlated to a different method of determining 
bone density. Thus, if bone density is measured using dual photon 
absorptometry, and one wanted to correlate the bone density coefficient to 
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10 



the results of the dual photon absorptometry, the overall relationship of the 
tests in the algorithm would be the same or similar as disclosed herein but 
the coefficients could be different 

In its preferred embodiment, the bone density probability quotient is 
assigned to one of the following diagnostic categories: 

Normal to mild osteopenia 
Moderate osteopenia 
Severe osteopenia 

In addition, the menopausal status may be determined based on the 
results of the estradiol and progesterone levels into one of the following 
diagnostic categories: 

15 Probable Pre-menopausal 

Probable Peri-menopausal 
Probable Post-menopausal 

Thus, according to the present invention, using blood concentrations 
20 of certain blood constituents, one can not only diagnose bone density but 

can also obtain an indication of the underlying cause of the osteopenia. 

This invention is further illustrated by the following examples, 
which are not to be construed in any way as imposing limitations upon the 
scope thereof. To the contrary, it is to be clearly understood that resort may 
25 be had to various other embodiments, modifications, and equivalents thereof 

which, after reading the description herein, may suggest themselves to those 
skilled in the art without departing from the spirit of the present invention 
and/or the scope of the appended claims. 

30 Example 8 

In the following example, the bone density coefficient algorithm was 
correlated with bone density as measured by radiographic absorptometry. 
The coefficients were calculated by correlating the blood concentration of the 
indicated blood constituents in 27 patients to the results of bone density in 
35 those patients as measured by radiographic absorptometry. The bone 

density algorithm which results is as follows: 
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0 + (-8.955)(Ca) + (79.370)(P) + (34.076)(E2) + (-9.216)(Pg) + 
(-0.600)(Alkp) + (-57.855)(I-Alkp) + (0.926)(CaCa) + (-2.735)(PP) + 
(-0.272)(E2E2) + (-0.064) (PgPg) + (0.004)(AlkpAlkp) + (-4.029)(CaP) + 
(-3.156)(CaE2) + (1.172)(CaPg) + (-0.031)(CaAlkp) + (8.238)(CaI-Alkp) 
+ (1.906)(PE2) + (-0.138)(PPg) + (-0.018)(PAlkp) + (-6.710)(PI-Alkp) + 
(0.193)(E2Pg) + (-0.048)(E2Alkp) + (-1.825)(E2I-Alkp) 

wherein: 

Ca = Serum concentration of calcium 

P = Serum concentration of phosphorus 

E2 = Serum concentration of estradiol 

PG = Serum concentration of progesterone 

Alkp = Serum concentration of total alkaline phosphatase 

I-Alkp = Serum concentration of serum lymphycyte-derived alkaline 

phosphatase isoenzyme. 

The bone density coefficient that is obtained from the bone density 
algorithm is assigned to one of three categories: 

Group 1 normal to mild osteopenia 

Group II moderate osteopenia 

Group III severe osteopenia 

Group IE £ approximately 78.5 < Group H < approximately 100 < Group I 
It has been determined that the bone density coefficient calculated according 
to the present invention places patients in the same groups as radiographic 
absorptometry measurements for bone density. Thus, the present invention 
provides a safe, economical and accurate method for diagnosing osteopenia. 

Bone density quotients derived practicing the present invention using 
the six biochemical serum constituents were found to strongly correlate with 
bone density measurements from radiographic absorptometry. It is to be 
understood that the preferred embodiment utilizes all six of the biochemical 
serum constituents listed hereinabove. However, if any five of the serum 
constituents are utilized and then included in the bone density algorithm, a 
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10 



15 



20 



30 



bone density quotient is obtained which does not correlate as well as when 
six tests are used. 

It is contemplated as part of the present invention that any five of the 
biochemical serum constituents can be used to determine the osteopenic state 
of a patient. For example, if serum calcium is removed from the bone 
density algorithm, then the algorithm predictive value for moderate or severe 
osteopenia is approximately 66% although the predictive value for normal to 
mild osteopenia is 100%. If progesterone is removed from the bone density 
algorithm, then the algorithm predictive value for moderate or severe 
osteopenia of approximately 75% and a predictive value for normal to mild 
osteopenia of approximately 94%. If alkaline phosphatase is removed from 
the bone density algorithm, then the algorithm predictive value for moderate 
or severe osteopenia of approximately 83% and a predictive value for 
normal to mild osteopenia of approximately 85%. If estradiol is removed 
from the bone density algorithm, then the algorithm predictive value for 
moderate or severe osteopenia of approximately 67% and a predictive value 
for normal to mild osteopenia of approximately 85%. However, it has been 
determined that when all six tests are included in the bone density algorithm, 
the predictive value for moderate and severe osteopenia is 100% and the 
predictive value for normal to mild osteopenia is 100%. All of the predictive 
values may vary slightly with different patient populations. 



25 woman: 



Example 9 

The following serum constituents are measured in a 37 year old 



Calcium 9.5 mg/dL 

Phos P hate 5.0 mg/dL 

Estradio1 2pg/mL 

Progesterone 0 .2 ng/mL 

Alkaline Phosphatase go u/L 

Alkaline phosphatase isoenzyme 0 (negative) 



35 



The above blood chemistry values were used to calculate a bone 
density coefficient using the following algorithm: 
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0 + (-8.955)(Ca) + (79.370)(P) + (34.076)(E2) + (-9.216)(Pg) + 
(-0.600)(Alkp) + (-57.855)(I-Alkp) + (0.926)(CaCa) + (-2.735)(PP) + 
(-0.272)(E2E2) + (-0.064)(PgPg) + (0.004)(AlkpAlkp) + (-4.029)(CaP) + 
(-3.156)(CaE2) + (1.172)(CaPg) + (-0.03 1 )(CaAlkp) + (8.238)(CaI-Alkp) 
5 + (1.906)(PE2) + (-0.138)(PPg) + (-0.018)(PAlkp) + (-6.710)(PI-Alkp) + 

(0. 193)(E2Pg) + (-0.048)(E2Alkp) + (-1.825)(E2I-Alkp) 
where: 

Ca = Serum concentration of calcium 
P = Serum concentration of phosphorus 
10 E2 = Serum concentration of estradiol 

PG = Serum concentration of progesterone 

Alkp = Serum concentration of total alkaline phosphatase 

I-Alkp = Serum concentration of serum lymphycyte-derived alkaline 

phosphatase isoenzyme. 

15 

A bone density coefficient of 101.2 is calculated for this patient. 
This coefficient falls within Group I in the following severity scale: 

Group in <, 78.5 < Group II S 100 < Group I 

20 

Group I normal to mild osteopenia 

Group II moderate osteopenia 

Group III severe osteopenia 

25 This patient has normal to mild osteopenia. When the patient's bone 

density is measured by radiographic absorptometry, she is found to have 
normal bone density. 
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Example 10 

The following serum constituents are measured in a 47 year old 
woman: 
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Calcium 9.5mg/dL 

Phosphate 5.0mg/dL 

Estradiol 2pg/mL 

Progesterone 0.2 ng/ml 

5 Alkaline Phosphatase 80 U/L 

Alkaline phosphatase isoenzyme l (positive) 



The values are inserted into the bone density algorithm recited in 
Example I and a bone density coefficient of 84.4 is calculated. This falls in 
Group II, indicating a moderate osteopenia. When the patient's bone 
density was measured by radiographic absorptometry, she was found to be 
moderately osteopenic. 

Example 11 

The following serum constituents are measured in a 47 year old 
woman who has been diagnosed with breast cancer: 



Calcium 9.0mg/dL 

Phosphate 4.0 mg/dL 

20 Estradiol Opg/mL 

Progesterone 0 ng/mL 

Alkaline Phosphatase 60 U/L 

Alkaline phosphatase isoenzyme 1 (positive) 



25 The values are inserted into the bone density algorithm recited in 

Example I and a bone density coefficient of 69.9 is calculated. This falls in 
Group m indicating a severe osteopenia insured by another method. 

Example 12 

30 ^ following example shows the correlation between the method of 

the present invention and Dual Photon Absorptometry (DPA) as a tool for 
diagnosing and determining the severity of osteoporosis in a particular 
individual. 200 female subjects representing a cross section of ages and 
menopausal status were evaluated. These specific age groups were 

35 represented as follows: 
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AgS Number of EatiSIttS 

25-35 26 

36-45 25 

46-55 33 

5 56-65 30 

66-75 34 

76-85 25 

85 and over 26 



10 Each subject completed an Osteoporosis Data Questionnaire Form, 

which included detailed information on personal history, family history, 
gynecologic history, medical history, surgical history, and drug history, 
with special references to hormone treatment. 

Three tubes of blood were collected from each subject via 

15 venipuncture. One tube was sent to Barnhill-MetPath Laboratories in 

Savannah, Georgia for analysis of serum calcium, phosphorous, and total 
alkaline phosphatase. The second tube was sent to Horus Therapeutics, 
Inc. in Savannah, Georgia for determination of lymphocyte alkaline 
phosphatase utilizing the Isoelectric Focusing Method. The third tube was 

20 s ^ n t to MetPath Laboratories in Teterboro, New Jersey for evaluation of 

estradiol and progesterone by immunoassay. Data from all three sources 
were then subjected to statistical analysis. At the time of the analysis of 
data, each subjects chart was critically reviewed for relevant clinical 
information. 

25 The Bone Mineral Density (BMD) of each subject was also 

measured by Dual-Photon Absorptometry using 153 Gadolinium as a source 
(using Lunar DP3, Lunar Radiation Corporation, Madison, Wisconsin). 
DP A measurements were performed on the Lumbar Spine (L1-L4) and the 
hip, specifically the femoral neck, Ward's Triangle, and the trochanteric 

30 regions. BMD was measured in gm/cm 2 for each region noted, and the 

fracture risk was determined. 

The BMD measurements were then age matched, and adjusted for 
sex, age, ethnic group, and weight. From the individual values obtained, 
average BMD for hip and spine was calculated. The overall average BMD 
35 for both hip and spine was also computed. 
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A mathematical algorithm correlating the results using the present 
invention with DP A was determined with the resultant data is shown in 
Figure 11. 

Figure 12 illustrates the discrepancies between the DP A 
5 measurements and the measurements obtained using the method of the 

present invention for individual subjects. The amount of the discrepancy 
(the residual) for any patient is the vertical distance between the points for 
that patient on the DPA and the present invention lines. 

Figure 13 illustrates the appropriateness of using the present 
10 invention as a prediction of DPA measurements. Note the absence of any 

tendency for points to be consistently above or consistently below the 
diagonal straight line, as greater DPA measurements are considered. 
Additionally, there is an absence of any tendency for points along any 
vertical line to be more widespread than are points along any other vertical 
15 line. Finally, note that the generally elongated shape formed by the points is 

consistent with a strong linear correlation between DPA and the present 
invention. 

Figure 14 demonstrates the normality of values predicted by the 
present invention (the horizontal scale). The vertical scale shows the of the 

20 normalized scale (z) scores values determined according to the present 

invention as would be expected from the assumption of the present 
invention values forming a normally distributed collection of numbers. 

The closer the dots are to the straight line, the greater is the 
consistency between what is observed about the present invention values 

25 and what is expected about them from the assumption of their forming a 

normally distributed collection of numbers. 

Figure 15 demonstrates the normality of residuals (the difference 
between the observed DPA measurement and the value predicated by present 
invention), as evidenced by the arrangements of points along the diagonal 

30 straight line. The horizontal scale shows values of residuals The vertical 

scale shows the residuals 1 normalized scale (z) scores as would be expected 
from the assumption that the residuals form a normally distributed collection 
of numbers. 

The closer the dots are to the straight line, the greater is the 
35 consistency between what is observed about the residuals and what is 
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expected about them from the assumption of their forming a normally 
distributed collection of numbers. 

Figure 16 demonstrates that the variance of residuals stays the same 
regardless of the value calculated according to the present invention involved 

5 in the residuals 1 computation. This homogeneity of variance appears in the 

uniform density of the dots in a rectangle formed by residuals between -0.17 
and +0.17 and by the present invention between +0.4 and +1.0. Moreover, 
there is no discemable tendency for points along any vertical line to be more 
widespread than are points along any other vertical line. 

10 Figure 17 demonstrates the independence of residuals and the values 

predicted by the present invention. The horizontal scales shows the values 
predicted by the present invention. The vertical scale shows the residual's 
studentized scale (t) scores. 

When there is no relationship between values predicted by the 

15 present invention and the residuals involving those predicted values (that is, 

when there is independence between residuals and predicted present 
invention values) the following criteria should be met: 

1 . There should be no discemable pattern (such as a straight or 
20 curved line) to the dots. 

2. About half the dots should appear in the top half of the 
graph. 

Figure 18 demonstrates that virtually all observed DPA 
25 measurements are satisfactorily predicted by the method of the present 

invention. The horizontal scale shows values predicted by the method of the 
present invention. The vertical scale shows Cook's distances between 

» 

observed DPA measurements and the predictions of those measurements by 
the method of the present invention Thus, the closer the dots are to a height 
30 of Cook = 0, the more accurate is the method of the present invention. 

Because Cook's distances tend to form a collection of numbers fitting an F- 
statistical distribution, the two outline points at heights between 10 and 20 
are not thought to be problematic in view of the large size of the sample. 
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Example 13 

The following example describes the training of a neural network to 
prognose prostate cancer 

As shown in Fig. 19, a total of 52 samples were divided into 2 

5 groups, a training set and a generalization testing set. The training set 

contained 40 samples (28 stable and 12 progressing) and the generalization 
testing set contained 12 samples (9 stable and 3 progressing). 

The initial network architecture was selected based on the level of 
complexity of the classification task. A multi-layer feedforward network 

10 was used. Selection of the initial architecture involved the selection of the 

number of hidden layers and the number of neurons in each hidden layer. 
Several trial iterations were performed to determine and adequate 
configuration that showed good results on both the training sample set and 
the generalization test sample set. The present network had one hidden 

15 layer, having nine neurons, and two output neurons. 

Initially, connection weights among the neurons were randomly set 
The neural network had five input neurons, corresponding to five input 
variables significant for prostate cancer: TPS, PSA, PAP, CEA, and 
testosterone. The training data is shown in Fig. 19. During training, the 

20 five input variables for each patient were first linearly scaled into the 

continuos range between 0.0 and 1.0. The resultant five numbers were then 
presented as an input vector to the input neurons of the artificial neural 
network. 

For each of the input vectors, the network generated an output based 
25 on the connection weights among the network neurons. The output can be a 

single value or a vector of numbers, depending on the number of output 
neurons used. The network used had two output neurons. The outputs of 
the two neurons were processed by the following mathematical equation to 
yield a single diagnostic index: 

30 

Index = (ANN2-ANN1) + 0.5 

2 

Each neuron in the network participated in the output calculation by passing 
35 the sum of all inputs to the neuron through a non-linear s-shaped function 

(often a ligistic function) and sending the result to each and everyone of the 
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neurons in the following adjacent layer. The generated output or each 
output neuron was compared to the desired "target" output. A value of 0.1 
corresponded to a diagnosis of stable and an output of 0.9 corresponded to a 
diagnosis of progressing. The difference was used to calculate an error term 
5 to guide the training algorithm, i.e., the back propagation algorithm, in the 

adjustment of network connection weights in an attempt to reduce the 
differences between network outputs and target values over the training 
sample set. 

After training, the neural network correctly classified 100% of the 

10 samples. 

When presented with the generalization test results, the trained 
neural network correctly identified 100% of the stable samples and 66% of 
the samples where the disease was progressing. 



15 It should be understood, of course, that the foregoing relates only to 

preferred embodiments of the present invention and that numerous 
modifications or alterations may be made therein without departing from the 
spirit and the scope of the invention as set forth in the appended claims. 
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Claims 



What is claimed is: 



1 . A method for diagnosing a disease in a human or animal 
5 comprising the steps of: 

measuring the concentrations of a predetermined set of biomarkers 
known to be associated with the disease from a biological fluid from 
the human or animal; 

10 

scaling the digitized values of the biomarkers concentrations; and 



introducing the scaled values to a trained neural network means, 
whereby the output values from the neural network means tend 
15 toward the upper value when the human or animal has the disease 

and the output values tend toward the lower value when the human 
or animal does not have the disease. 



2 . An apparatus for diagnosing a disease in a human or animal 
20 comprising: 

means for digitizing the concentrations of a predetermined set of 
biomarkers known to be associated with the disease from a 
biological fluid from the human or animal; 

25 

a means for scaling the digitized values; and 



a trained neural network means coupled to the digitizing and scaling 
means for generating network output values between an upper and 
30 lower value, whereby the output values tend toward the upper value 

when the human or animal has the disease and the output values tend 
toward the lower value when the human or animal does not have the 
disease. 
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3. A method for diagnosing a disease in a human or animal 
comprising the steps of: 

measuring the concentrations of a predetermined set of biomarkers 
known to be associated with the disease from a biological fluid from 
the human or animal; 

scaling the digitized values of the biomarkers concentrations; 

introducing the scaled values to a first trained neural network means, 
whereby the output values from the first neural network means tend 
toward the upper value when the human or animal has the disease 
and the output values tend toward the lower value when the human 
or animal does not have the disease; and 

introducing the output value from the first neural network and a 
second set of predetermined biomarkers to a second trained neural 
network means, whereby the output values from the second neural 
network means tend toward the upper value when the human or 
animal has the disease and the output values from the second neural 
network means tend toward the lower value when the human or 
animal does not have the disease. 

4. The method of Claim 3 wherein the second set of 
predetermined biomarkers includes patient biographical information. 

5. The method of Claim 3 wherein the second set of 
predetermined biomarkers are the same as the first set of predetermined 
biomarkers. 
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0.825 


549 


26 


9.5 


3.9 


50 


0.4 


46 


7.2 


49 


0.862 


560 


33 


8.9 


3.4 


118 


0.3 


79 


7.4 


61 


0.791 


562 


23 


9.4 


2.2 


-1 


0.4 


58 


0 


36 


0.815 


579 


20 


9 


3,4 


48 


0.5 


90 


0 


53 


0.797 


588 


36 


9.3 


3.5 


-1 


0.1 


48 


0 


51 


0.793 


591 


2? 


9.7 


3,7 


36 


0,5 


52 


0 


54 


0.93? 


593 




9 


3.7 


108 


0.3 


69 


9.6 


41 


0.81 1 


KOI 






■?TY 








0 


71 




602 


49 


9.3 


3.6 


180 


4.3 


69 


0 


46 


0.693 


2520 


29 


9.2 


4.4 


» 


0.2 


49 


0 


41 


0.877 


3000 


28 


9.4 


3.7 


49 


0.4 


66 


0 


53 


0.828 


3020 


25 


9.6 


4 


164 


16.8 


56 


0 


41 


9880 


100 


52 


8.7 


3,1 


?7J 


0.3 


59 


o 


54 


0.703 


113 


53 


9.8 


3.5 


3 


0,? 


63 


0 


67 


0.690 


126 


42 


9 


4 


99 


16.4 


8? 


14.4 


36 


0.984 


130 


50 


9.3 


3.5 


336 


0.6 


40 


0 


25 


0.754 


139 


74 


8.9 


2,7 


209 


0.4 


35 


5,9 


36 


0 648 


160 


60 


9.5 


3.9 


81 


0.3 


47 


5.2 


53 


0 "04 


176 


_56 


9 


2.9 


115 


0.2 


64 


o 





0 677 



m =7 

SUBSTITUTE SHEET (RULE 26) 



WO 96/12187 



11/24 



PCT/US95/01379 



| PATIENT 


AGE 


CAL 


PHQS 


ETWQ 


PROG 


TALP 




PLIVER 


SCORE 


H 177 


39 


9.4 


3.1 


71 


0.2 


82 


0 


45 


0.751 


208 


42 


9.1 


2.4 


177 


15.6 


54 


0 


70 


0.757 


II 211 


50 


9.5 


3.1 


69 


0.2 


77 


0 


60 


0.696 


II 214 


25 


9,2 


3.7 


-1 


Q.9 


4? 


0 


29 


0.852 


H 251 


24 


9.5 


3.9 


-1 


0.1 


61 


0 


25 


0.844 


H 273 


3Q 


9,4 


3 


99 


0.2 


156 


8.3 


26 


0.776 


II 299 


27 


9,2 


3,2 


41 


0.2 


77 


0 


58 


0.788 


|| 309 


35 




3.9 


55 


0.3 


90 




62 


0.773 


H 316 


41 


9.4 


3.5 


75 


7.3 


44 


0 


27 


0.733 


|| 320 


52 


9.3 


3.4 


57 


0.2 


58 


0 


55 


0.707 


II ? 55 


54 


9.1 


3.6 


93 


0,4 


55 


0 


40 


0.718 


II 358 


54 


9.5 


3,6 


38 


0.3 


71 


0 


54 


0 688 


II 361 


49 


9 


3.8 


31 


0,3 


74 


9.5 


50 


0.744 


II 37 1 


49 


9.5 


3.8 


-1 


0,2 


7? 


10.3 


42 


0.698 


II 377 


52 


9.6 


3.6 


120 


0.2 


68 


0 


5? 


0.717 


II 385 


55 


9 


3,3 


-1 


0,2 


97 


0 


73 


0.646 


|| 398 


53 


9.3 


4.6 


30 


0.2 


97 


10.3 


28 


0.697 


|| 406 


26 


9.6 


3.4 


35 


0.3 


62 


0 


70 


0.806 


II 414 


52 


9,3 


3.9 


56 


0.4 


66 


0 


72 


0.710 


421 


31 


9 


3.5 


72 


0.3 


69 


0 


52 


0.818 


536 


48 


9,4 


3,4 


54 


3.2 


56 


0 


63 


0.688 


1 544 


38 


9.3 


2.4 


126 


0.2 


64 


0 


59 


0.747 


S66 


36 


9,3 


4 


-1 


0.5 


41 


0 


44 


0.822 


II 603 


27 


8.9 


4.4 


42 


0.3 


48 


0 


38 


0.872 


II 109 


54 


9.4 


3.3 


64 


0.3 


121 


9.6 


65 


0.652 


1 110 


58 


9.4 


3.9 


20 


0.1 


85 


0 


58 


0.663 


II 129 


45 


8,8 


3.7 


197 


10,8 


56 


0 


53 


0.733 


1 141 


46 


9.7 


3.6 


141 


2.1 


68 


0 


44 


0.709 


1 158 


45 


8.9 


3.3 


243 


0.8 


63 


0 


43 


0.763 


1 164 


51 


9,3 


3.5 


139 


0.5 


59 


0 


43 


0.727 


1 170 


51 


9.4 


3.3 


120 


4.2 


50 


0 


30 


0.737 


1 178 


36 


9.7 


3.3 


139 


0.4 


40 


0 


46 


0.824 


II 180 


38 


8,7 


3,4 


48 


0.4 


52 


0 


39 


0.783 


II 181 


34 


8.7 


3.4 


78 


U.J 


67 


6.8 


33 


0.817 


1 183 


44 


9 


3.6 


131 


0.2 


35 


0 


54 


0.807 


1 191 


45 


9.3 






fl A 


32.., 


o 






ll 194 


29 


9.2 


3.7 


89 


0.1 


82 


0 


52 


0.824 


1 vrn 


45 


9.2 


3.7 


226 


1.6 


46 


o 


27 


0,786 


204 


25 


9.5 


3.5 


170 


0,3 


4? 


o 


46 


0-863 


209 


42 


9.6 


2.9 


539 


0.4 


50 


0 


64 


0.738 


216 


46 


9.2 


3.6 


156 


0.5 


43 


0 


39 


0.780 


218 


43 


9.5 


4.4 


133 


0,7 


59 


0 


41 


0.791 


225 


37 


9 


2.7 


159 


1.4 


64 


6.7 


48 


0.766 


226 


72 


9.9 


4.9 


46 


0.3 


77 


8.9 


61 


0.643 


227 


46 


9.6 


4.4 


-1 


9,2 


94 


0 


38 


0.757 




43 


9,2 


3,4 


94 


0 


71 


0 


77 


0.721 




2£- 




3-2 


-1 


0.3 


51 


0 


61 


0.777 




SUBSTITUTE SHEET (RULE 26) 



WO 96/12187 



12/2 4 



PCT/US95/01379 



K 



PATIENT 




CAL 


r»ur/~\c 






TALr 


TTXITWC 


or rvrpp 




711 


5fi 


Q 1 


1 5 


41 


0 3 


73 


^^^^ 


23 


0 709 


7ifl 
Z38 


4A 
*ry 


52 i 

O.J 


1 4 


Old 


0 7 




5 0 


46 ! 


1 094 


7 in 


4V 


Q 5 


1 ft 


41 


1 ft 


71 


0 


63 


0 758 


7/iz; 


2±_ 


O 52 


1 1 


A 


0 ^ 
Y ir 


47 


n 


58 

JO 


0 804 


74Q 


31 


52 0 


1 52 


7ft 


0 1 


ftl 


o 


62 


0 724 




1 O 


O A 


4 1 


40 


0 ^ 

^^^^ 


52ft 




34 


0 841 


7<7 


71 


52 O 


7 0 


1 152 


7 7 


75 


o 


21 


0 816 


7A7 


7A 


O 1 


1 ft 

a us 


.1 
f 


0 7 

^^^^^ ! 


47 


0 1 


24 


0 844 


7ft7 


71 


Q 


^ ft 


170 ! 


0 4 


84 ^ 


o 


50 


0.825 


77fi 


7 1 


Q 1 


3 9 


47 


0 9 


40 


o 


70 ! 


0.850 




77 ! 


0 4 


4 ft 


58 


0 2 


138 


o 


69 


0 811 


7751 


70 


0 4 


3 Q 


404 


1 2 


62 


7 9 


46 i 


0.864 


7C7 


7ft 


Q 4 


1 0 


ftO 


0 2 


57 


o 

^^^^ 


36 


0 859 


7S1 


75 




^ ft 


_1 


0 3 
^^^^^^ 


48 


o 


32 


0 847 


75252 


74 


0 7 


4 1 


65 


0 3 


107 


7 2 

IT* 


35 


0.816 


7 CO 


77 
s L 


0 4 


3 4 


68 


0 6 


87 


0 i 

^^^^ 


37 


0.812 


701 


7ft ! 


0 ft 


3 3 ! 


.1 
1 


0 3 


58 


o i 


28 i 


0.823 


7Q1 


7ft 


Q A 


4 


32 


0 2 


68 


0 

^^^^ 


35 ! 


0.838 


705 


7ft 


Q 1 


3 3 


57 i 


0 5 


51 


o 

^^^^^ 


63 


0.324 


n n. 


47 


O 5 


^ ft 




1 ^ 


53 


o 


46 


0 730 


114 
314 


452 


O 1 

9,1 


7 7 
it* 1 


_i 
i 


0 2 


64 


o 


53 S 


0 704 


77fl 


17 

3/ 


0 4 


1 7 


54 


0 4 

V .*t 


74 


o 


55 


0 762 


IK 


47 


O 7 

9.2 


4 7 


40 
ZZ 


0 4 


74 


o 


56 


0 745 


1A 1 


17 


O 1 
7iP 


1 1 


.1 

-L 


0 1 


ftft 


o 


56 


0 758 


^ A ft 

349 


50 


Q Q 


1 1 

3.1 


151 


n t 

Yt * 


ft7 

D / 


52 
o 


30 


0 710 


o f i 
353 


5Q 


9 


A 1 

4.1 


i7vi 


n 7 


f ^1 


A 

u 


i io 

iLL 


0 750 


. 354 


58 


Q 7 

9.7 


4 1 


74 


n i 


452 


ft 8 

Vj . O 


! 45 


0 726 


1 222- 




O 1 


1 1 

3-3 


7ft 


0 ft 


01 


13 8 


31 


0 694 


3oo 


7>4 

Z4 


O 4 

9.4 


7 7 
2./ 


.1 


o s 


50 


o 


31 


0 806 


372 


7 1 


1 A 


7 O 

2.y 


5A 


O 7 


70 


— — LL_ 


61 


0 789 


I 3/5 


A1 

42 


O 7 


1 ^ 

3-5 


1ftO 


n i 


77 


n 


56 


0 754 


1 ioi 


62 


9-5 


1 A 




n i 


ftl 


0 

- v 


30 


0 672 


?9 4 


A7 


y.4_ 


4 


14 




0 7 


08 


! 75 3 


33 


0 644 


1 ?9 ' 


<7 




1 7 


152 


n 7 


111 




55 


0 669 


1 388 


64 


9.1 


3.1 


87 


0.3 


80 


0 


75 


0.636 


1 396 


57 


8.8 


3.1 


76 


Qi? 


83 


7,7 


51 


0.659 


1 397 


58 


9,6 




73 


0.3 


85 


0 


52 


0.665 


400 


50 


9.1 


3,5 


1 


9,2 


99 


9.7 


50 


0.671 


410 


37 


9,5 


2.8 


-1 


0.3 


44 


0 


48 


0.769 


1 411 


49 


9.4 


3.1 


138 


7.5 


69 


o 


27 


0.674 


1 41? 


31 


9.6 


4.6 


-1 


0.4 


42 


0 


72 


0.861 


424 


54 


?,4 


3.1 


59 


0.4 


101 


Q 


35 


0.662 


1 597 


45 


9 


3.4 


46 


0.3 


55 


o 


65 


0.745 


I 513 


25 


9,4 


3.5 


52 


0.2 


72 


0 


52 


0.815 


1 533 


61 


9,6 


3.9 


-1 


0.1 


80 


0 


61 


0.655 


1 535 


24 


9 




-1 


0.3 


68 


0 


55 


0.795 


1 537 


48 


9.9 


3.2 


32 


0.3 




0 


58 


0.679 
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SUBSTITUTE SHEET (RULE 26) 



WO 96/12187 



U /2 4 



PCT/US95/01379 



PATIENT 


AGE 


CAL 


PHOS 


ETWO 


PROG 


TALP 




PLIVER 


SCORE 


317 


47 


9.6 


2,4 


-1 


0.2 


60 


0 


27 


0.709 


324 


46 


9.4 


4 


117 


P,4 


75 


0 


54 


0.732 


327 


46 


9 


2,4 


322 


1 


46 


0 


52 


0.765 


329 


49 


9.5 


3,7 


72 


2.4 


78 


0 


45 


0.692 


337 


46 


9.3 


3,5 


122 


19,2 


48 


0 


42 


0.700 


352 


56 


8,8 


3.6 


73 


0,2 


53 


0 


35 


0.707 


357 


35 


9,4 


2.9 


89 


2,2 


71 


7.6 


23 


0.791 


367 


42 


9.1 


3.2 


58 


0.5 


90 


9.8 


64 


0.846 


369 


62 


9.6 


3.6 


31 


0.2 


86 


0 


39 


0.645 


379 


42 


8.6 


3.4 


156 


6,2 


84 


0 


61 


0.742 


380 


49 


9.6 


4.2 


-1 


0.2 


71 


11 


52 


0.690 


383 


62 


9 


2.5 


81 


0,1 


66 


7 


52 


0.665 


392 


63 


9.5 


3.4 


-1 


0.3 


98 


0 


76 


0.620 


517 


25 


10 


4 


31 


0.1 


50 


0 


50 


0.853 


526 


47 


9.1 


3.8 


94 


9,2 


« 


0 


51 


0.750 


543 


59 


?,4 


4.2 


-1 


0.2 


104 


o 


20 


Qt$62 


553 


48 


8.8 


4.1 


35 


0,5 


57 


0 


66 


0.759 


556 


69 


9.6 


4 


77 


0.1 


44 


6.8 


53 


0.668 


567 


29 


9.4 


4.2 


-1 


0.3 


67 


9.7 


2? 


0.841 


575 


45 


9,1 


3.1 


406 


0.6 


104 


0 


54 1 


0.704 


584 


57 


9.2 


3.2 


141 


0.2 


79 


0 


40 


Q T $74 


599 


49 


8.9 


4 


42 


0.2 


79 


9 


51 


0.748 


2510 


26 


9.5 


3.9 


157 


16.2 


56 


0 


25 


0.881 


105 


75 


9.4 


4.1 


42 


0,4 


88 


8 


52 


0.605 


II 116 


67 


8.9 


4 


-1 


0.3 


56 


0 


70 


Q.$47 


II 118 


55 


8.9 


4.3 


38 


0,? 


78 


0 


37 


0.699 


II 119 


56 


9,3 


4.3 


75 


0.3 


41 


0 


43 


0.740 


127 


73 


9 


4.1 


30 


0.4 


54 


6.5 


35 


0.632 


134 


50 


9.4 


3.7 


486 


0.4 


54 


0 


6? 


0.737 


136 


36 


9,1 


3.6 


167 


5.7 


154 


0 


74 


0.758 


140 


56 


9 


3.9 


-1 


0.2 


94 


0 


58 


0.663 


143 


55 


9.1 


4 


-1 


0.2 


65 


0 


38 


0.698 


145 


49 




3.3 


67 


8.9 


63 


0 


54 


0.656 


II 147 


53 


9.1 


3.6 


75 


0.2 


103 


4.9 


78 


0.670 


II 157 


52 


9.2 


3.5 


65 


9.2 


70 


0 


62 


! 0.701 




en 

5 ' 


O A 


•* 


i 




1 7Q 


30.2 


i 77 




161 


74 


9.3 


3.9 


-1 


0.1 


62 


o 


34 


0.623 


162 


62 


9 


4.1 


81 


0,1 


54 


7.2 


21 


0.698 


163 


46 


9,2 


2.9 


128 


5,7 


62 


0 


35 


0.704 


165 


57 


9,? 


2.9 


-1 


u- 


108 


0 


39 


0.629 


169 


59 


9.5 


3,9 


61 


0,1 


64 


0 


34 


0.693 


174 


40 


9.5 


3.2 


261 


4.6 


85 


0 


2? 


0.764 


175 


39 


9.1 


2.3 


126 


3.5 


90 


0 


36 


0-726 


185 


78 


8.7 


2.9 


-1 


0.1 


99 


8.5 


35 


0.563 


186 


73 


9,2 


2.9 


62 


0.2 


55 


0 


42 


0.616 


187 


59 


9 


3.8 


55 


1.8 


70 


0 


37 


0.682 


190 


41 








116 




r^n °- 


2$_ 


0.772 




io 



SUBSTITUTE SHEET (RULE 26) 



WO 96/12187 



15/2/, 



PCT/US95/01379 



PATIENT 
192 
195 
202 
210 
215 
221 
224 
229 
237 
241 
263 
294 
308 
313 
319 
322 
323 
334 
338 
339 
343 
344 
351 
362 
373 
376 
390 
395 
399 
401 
413 
423 

527 

552 

551 

557 
563 
565 
570 
573 

521 

577 
585 
590 
594 

601 

601 



AGE 

42 

-JJL 

67 

_£0_ 
_45_ 

52 

_52_ 

_5J_ 
48 

_2i_ 
_25_ 

_54_ 

53 
31 
42 
56 

_6S_ 

42 

_52_ 

44 

61 
46 
53 

_62_ 

70 

_22_ 
_51 

49 
83 
46 
44 
38 

_54 
50 
59 
44 

__7p_ 
24 
65 

_46_ 

_62- 

65 

^55. 

45 



CAL 

_L5. 
__2_ 
_24_ 

9.6 

_2i 
_L5. 

8.7 
9.4 

JSJl 

8.8 

9. 

9.4 
8.8 
8.8 
9.9 
9.9 
9.1 
9.2 
9.2 
9.6 
9.9 
9.9 
9_ 

-Si 

9.4 
8.7 
_L4 

9.2 

-Li 

9.4 

_£L5. 

9.4 
9.5 



9.2 

£2. 

9.2 

9 

22. 

8.9 
9.7 



PHOS 
4.1 
3-4 
3.2 

3- 2 

3,8 
4_ 

3.6 

4- 1 
4_ 

3.9 

4.2 

3- 8 

3,1 
4.3 

3.5 

3.3 

3.6 

4.1 

3.2 

3.7 

3.4 

4.3 

3.9 

3.9 

3.4 

4- 7 
3.8 
4.1 

1 

3.8 
3.6 
3-6 
3.3 

3- 7 
3_ 

3.5 
3.8 
3.2 
3.5 

4,2 

3.2 

4- 1 

3,? 
3-7 
3.6 

3-1 



ETWO 

4£_ 

236 
66_ 

S2_ 

82_ 

no 

22. 

-\_ 

44 

22. 

2L 

-J_ 

j_ 

296 

-J_ 

2S_ 

S3_ 

25. 

84. 

55. 

82. 

140 

21. 

172 

iL 

J. 

iL 

-_L 

21 

20. 

62- 

40. 

d. 

108 
192 

21 

-± 

50. 

-± 

122 
397 

56. 

zL 

J. 

Si 

-l 



PROG 
04 

0.1 
0.3 

01 

1- 8 
0.8 

0-2 
0.3 

0-2 
0.5 

0-5 

2- 9 

0-2 
0.4 

2- 8 

0-3 
0.1 

9,4 

L 

0- 2 

1- 3 
0.3 

0-3 

0-2 
0.1 
0.2 

0-3 
0.4 

0.2 

9,1 

0-2 
0.4 

0-2 
0.2 

0-2 
0.5 
0.2 
0.5 

03 
0.3 

.04 

0.2 
0-2 
02 

3- 4 
0.1 



TALP 

5S_ 

61 

59 
49 

70_ 

60 

52 

97. 

90 

68 
69 

64. 

40_ 

22. 

49_ 

63 

54. 

44. 

22. 

78_ 

75 
116 
58 
83 

SI 

109 
108 
96 
64. 

HI 

22. 

47 
62 

SP_ 

126 

42 

46_ 

54. 

53 

52. 

64. 

26. 

77 
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* * Program 
** Author 

* * Date 

* * Purpose 



Calc. prg 

Terminal solutions (JJK) 
14 August 1992 
Quios calculating program 



set decimal to 16 
set fixed on 

****** Set Coefficient information to valuables ****** 



store 
store 
store 
store 
store 
store 
store 
store 
store 
store 
store 
store 
store 
store 
store 
store 
store 
store 
store 
store 
store 
store 
store 
store 
store 



0.000000000 to outnumb 
1009.33 to constant 



0.014786248 

424.462000000 

66.3468 

0.025761569 

0.019715699 

0.027232306 

0.000000263 

12.0319 

7.499334732 

2.220366330 

0.313008013 

0.016909117 

0.001369629 

0.000428831 

0.001798351 

0.001841052 

0.007649946 

0.002153604 

4.606320653 

0.119876531 

0.011786744 

0.000170224 

0.000000706 

outnumb= (constant) 

outnumb=outnumb + 

outnumb=outnumb + 

outnumb=outnumb + 

outnumb=outnumb + 

outnumb=outnumb 

outnumb^ outnumb 

outnumb-outnumb 

outnumb-outnumb 

outnumb=outnumb 

outnumb-outnumb 

outnumb-outnumb 

outnumb-outnumb 

outnumb=outnumb 

outnumb=outnumb 

outnumb-outnumb 

outnumb-outnumb 

outnumb=outnumb 

outnumb-outnumb 

outnumb=outnumb 

outnumb-outnumb 
outnumb-outnumb 

outnumb=outnumb 

outnumb=outnumb 



+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 
+ 

+ 

+ 




ML 



to c _age 
to c_ cal 
to c_cal2 
to c_csage 
to c_csphos2 
to c_cstintes 
to c_etwoimx2 
to c_ phos 
to c_ phos2 
to c_ phos3 
to c_ phos4 
to c_ phos5 
to c_ prog2 
to c_ progpliv 
to c_ progtint 
to c_ tnetwoimx 
to c_ tntintes 
to c_ agephos 
to c_ ca!3 
to c_ ca!4 
to c_ talp 
to c_talp2 
to c_talp3 

(c_age * age) 
(c_ cal * cal) 
(c_ cal2 * (cal * cal)) 
(c_ csage * (cos(age))) 
(c_ csphos2 * (cos(phos * phos))) 
(c_ cstintes * (cos(tintes))) 
(c_ etwoimx2 * (etwoimx * etwoimx)) 
(c_ phos * phos) 
(c_ phos2 * (phos 
(c_ phos3 * (phos 



* 

* 



* phos)) 
* phos. 



(c_ phos4 
(c_ phos5 
(c_ prog2 
(c_ progpliv * 
(c_ progtint 4 
(c_ tnetwoimx 
(c_ tntintes * 



* phos 

* phos 



phos)) 



(phos * phos 
(phos * phos 
(prog * prog)) 

* (prog * pliv)) 

* (prog * times)) 
(tan (etwoimx))) 
(tan(tintes))) 



* 
* 



phos)) 
phos * 



phos)) 



(c_ agephos * (age * phos)) 

(c_ cal3 * (cal * cal * cal )) 

(c_ cal4 * (cal * cal * cal * 
(c_ talp * talp) 

(c_taip2 * (talp * talp)) 

(c_ talp3 * (talp * talp * talp)) 



cal)) 
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N 



i = 0_N-1 4 input nodes 



Nl 8 



j = O NI 1 8 neurons in first hidden layer 



N2 1 



1 output neurons 



f(x) 



1 



1+e 



Logistic neuron output function 





100 




0 


I max. = 


25 


lmin. = 


10 




40 




0 




40 




0 



I range. = Imax. -lmin 



Wl = READPRN (link l_wgt) 



W2 = READPRN (link 2_wgt) 



Wl = 
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9x(5+l) weitht matrix between input 
layer and the first hidden layer 



W2 = (0.43 1.46 1.39 3.91 3.96 2.16 2.25 -6.36 -4.1) 2x9 weight matrix between 

hidden layer and output 
layer 
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